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1.0 INTRODUCTION 

During the final year of Grant sponsorship a paper reporting work 
done previously was published (Kidder, et al., 1980); a paper covering 
part of the work done this year was published (Chan et ai., 1980), and 
a Master's thesis was completed (Durkee, 1980). All of these papers 
are reproduced in Appendix B. A complete list of puclications sponsored 
by Grant NSG-5258 is presc.,«:ed in Appendix A. Scientific results ob- 
tained during the period 1 April 1979 to 31 March 1980 are summarised in 
the next section. 

2.0 SUMMARY OF SCII-NTIFIC RESULTS 

Work during this period concentrated on three areas: forecasting 

intensity change in tropical cyclones, forecasting tropical cyclone 
turning motion, and detecting summer precipitation over the North 
Central U.S. 

2.1 Tropical Cyclone Intensity Change 

It ha'j often been mentioned that for tropical cyclones to form or 
intensify, they must exist in large-scale environments of low-level posi- 
tive relative vorticity and upper- level negative relative vorticity 
(Riehl, 1948, 1950; Yanai, 1961; Dvorak, 1975; Sadler, 1976, 1978). 
Recently McBride (1979) has done compositing studies of developing versus 
non-developing c'.oud clusters. He found that the single most important 
indicator of tropical cyclone genesis was the difference between low- 
level (900 mb) and upper-level (200 mb) relative vorticity averaged over 
the area 0°- 6° latitude centered on the precursor cloud cluster. 

*^^900 " ^200^ 0°-6° larger by a factor of three for developing 


systems. 
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By thermal wind arguments, this vorticity difference is related to 
the laplacian of the mean temperature of the 900 mb to 200 mb layer: 


^900 


’200 


f '•200-' 900-200 


Cl) 


where R is the gas constant, and f the Coriolis parameter,. When averaged 
over the 0°-6° area, this becomes 

^900-200 0°-6° ■ ~ f^^^OOl r,o *-3r ^900-200^.0 

C) o 

where square brackets indicate averaging around the circle wlicse radius 
is indicated by the subscript, and r^o is the radius of the 6° latitude 
circle. Eq. (2) states that if a tropical disturbance is going to de- 
velop (intensify) one should see an inward gradient of I’goo-'^OO' 

To see whether satellite sounder data may play a role in forecasting 
tropical cyclone intensity, radial profiles of mean 1000 mb to 250 mb 
temperature from the Nimbus 6 Scanning Microwave Spectrometer (SCAMS) 
were constructed around eight intensifying tropical stonns in the western 
Pacific during the period August, 1975, to January, 1976 (Fig« 1). Seven 
showed distinct im^ard temperature gradients, required for intensification. 
The eighth (Alice) did not display an inward gradient and indeed was 
decaying 24 h later. 

Because the current operational polar-orbiting satellites (TIROS-N 
and NCAA-6) provide a dense network of soundings in the tropics, the 
above results indicate that they may play a key role in forecasting in- 
tensity changes in tropical cyclones. 

2.2 Tropical Cyclone Track Forecasting 

In order to better understand the relation between a tropical cy- 
clone's motion and the mass and flow fields of its environment, J. C. . L. 
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Chan and W. M. Gray of Colorado State University composited rawinsonde 
observations around West Indies tropical cyclones for the period 1961- 
1977. Wlien they stratified their data into three classes — storms which 
turned right, storms which turned left, and storms which continued on 
a straight track — they found that 24 to 36 h before a storm turned, the 
vertical shear C<ii^ference betva /a 300 and 200 mb) of the component of the 
environmental wind in the direction of storm motion (V^) was distinctly 
different in the three cases. For storms which turned right 24 to 36 h 
later, Vp decreased with height; for storms which turned left Vp in- 
creased wi.th height; and for straight-moving storms Vp changed little with 
height. Although the exact dynamical connection between vertical shear 
of Vp and storm direction is not fully understood, it appears to have 
potential as a forecast tool. 

By thermal wind considerations the shear of Vp should be reflected 
in the horizontal gradient of the mean temperature between 900 and 200 
mb CTgQQ_ 2 QQ) such that turning storms turn toward cold air. Because 
microwave sounders can estimate such mean temperatures even in the pres- 
ence of clouds, we were asked by Chan and Gray to investigate the pos- 
sibility that satellite data might be used to forecast tropical cyclone 
turning motion. Estimates of the mean 1000 to 250 mb temperature field 
around eleven tropical storms in 1975 were obtained from Nimbus 6 SCAMS 
data tapes. Analysis of these data showed that for turning storms, in 
all but one case, the turn was signaled 24 h in advance by a significant 
temperature gradient perpendicular to the storm's path at a distance of 
9° to 13° in front of the storm. Similarly, in all but one case, storms 
which continued on a straight course showed no significant gi'adient of 
T'iooo 250“ results have been submitted for publication (Chan et al.. 
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1980; see Appendix B) . They may have the potential to improve tropical 
cyclone 24 h position forecasts, the error in which is cui'rently nearly 
200 km. 

2.3 Summer Precipitation Frequency 

The lUectrically Scanning Microwave Radiometer on board Nimbus b 
(ESMR-6) measures both vertically and horizontally polarized upivolling 
radiation in a 250 MHz band centered at 37 Gflz. At those frequencies, 
precipitation appears cold in comparison with high omittance, high 
temperature land. A thresholding technique was developed to distinguish 
cold i-ain from warmer, non-raining areas. Rain was separated from 
standing water by examination of the two polarizations. Fig. 2 shows 
a comparison betw ' 'n ESMR-6-derived area of precipitation and raingage 
measurements . The thresholding technique was apiJlied to the North 
Central U.S. during the Summer of 1976 to estimate jirecipitation fre- 
quency. Good agreement was found between gauge measurements and satel- 
lite estimates except in Minnesota where a large number of small lakes 
probably biased the estimates. 

3.0 SUGGESTED FUTURE WORK 

One half of the earth's surface lies between 30°N and 30°S, and 
47 percent of that area is covered by ocean where few conventional 
meteorological observations are made. Satellite observations, tliere- 
forc, liave the potential to substantially improve weather analysis and 
forecasting in the tropics . 

Satellite sounder data has not been used extensively in the tropics 
because the uncertainty in the retrieved temperature at a particular 
level is often greater than the climatological variation. During the 
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Figure 2. Satellite derived rain areas and amounts of rainfall 
(hundredths of an inch) at surface observing rain- 
gages for 14 June 1976 at 0630 GMT. Area of satel- 
lite observation is also indicated. 
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course of tliis Grant, however, wo have shown that if one carefully 
chooses the satellite data and the meteorological situation to be ana- 
lyzed, satellite data can be very beneficial, IVe strongly recommend 
that efforts to apply satellite data to tropical weather analysis and 
forecasting be increased. 

Significant temperature gradients occur over the tropical oceans 
only near tropical cyclones. During the first year of Grant sponsor- 
ship we showed that brightness temperature anomalies in the 55.45 GHz 
cliannel of tlie Nimbus fa SC.\MS (caused by upper tropospheric temperatuer 
anomalies) could be used to estimate surface wind speeds outside of the 
radius of maximum winds. Because this technique is suitable for use 
with the TIROS-N Microwave Sounding Unit (M3U) or with tlie Defense 1 

Military Satellite Program (DMSP) microwave sounder (SSM/T) , we believe 
the technique should be tried with those data. 

During the second year of Grant sponsorsnip, work concentrated on 
using mean tropospheric temperatures to forecast 1) turning motion of 
tropical cyclones, and 2) intensity change. Initial results of these 
investigations (described above) look extremely promising. We highly 
recommend that investigation in these areas continue, and one of us 
(S. Q. K.) has proposed to do so. 
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ABSTRACT 

Upper tropospheric temperature anomalies arc detected in brightness temperature data from the 
Nimbus 6 Scanning Microwave Spectrometer (SCAMS). Brightness temperature anomalies are related 
to surface pressure anomalies through the radiative transfer and hydrostatic equation. Surface wind 
speeds at outer radii are then estimated using the gradient wind equation and a shearing parameter. The 
method is lirst tested using simulated satellite data constructed from temperature, pressure and height 
data recorded by aircraft reconnaissance of four hurricanes. Wind speeds in the 80^95 kl^a region arc 
estimated with 2-3 m s ‘ accuracy. Next, 55.45 GHz SCAMS data over eight typhoons during 1975 are 
used to estimate the radii of 15.4 m s ‘ (30 kt) and 27.5 m s ' (50 kt) winds. Accuracies of about +80 and 
+70 km, respectively, are found. It is suggested that the technique be further tested using data from the 
Microwave Sounding Unit (MSU) on board the TIROS-N and NOAA 6 satellites. 


1. introduction 

Over large areas of the oceans, the only feasible 
method for monitoring tropical cyclones is through 
the use of satelliteborne instrumentation. Since the 
launch of the first meteorological satellite in 1960, 
scientists have been collecting such data, Today the 
data are abundant and the challenge is to translate 
the raw data into accurate estimates of meteorologi- 
cal parameters. 

One of the most important parameters of a tropi- 
cal cyclone is surface wind speed. Among the uses 
of surface wind speed information are the following. 
Numerical models of storm motion need the wind 
speeds as input. Those responsible for routing ships 
need to chart courses which are outside the radius 
of 30 kt winds. Storm surge models are very sensi- 
tive to wind speed. Over the years a great deal of 
work has gone into the estimation of maximum 
sustained surface wind speed (intensity). The cul- 
mination of these efforts is the widely used Dvorak 
(1975) technique which is based on a correlation of 
intensity with cloud patterns. Unfortunately, the 
Dvorak technique does not estimate wind speeds 
at outer radii which are poorly related to intensity 
(see Gray and Frank, 1978, Tables 14 and 15). Re- 
cently, however, because of their importance to 
ship routing and storm surge prediction (Jelesnian- 
ski and Taylor, 1973), outer winds have become a 
topic of interest. Rodgers et al. (1979) have obtained 
good estimates of low-level wind speeds around 


‘ Current affiliation: Laboratory for Atmospheric Research, 
University of Illinois, Urbana 61801. 


tropical cyclones by tracking clouds in high-fre- 
quency (7,5 min) GOES (Geosynchronous Opera- 
tional Environmental Satellite) images. Initial re- 
sults from the Seasat-A Scatterometer System 
(SASS) and the Scanning Multichannel Microwave 
Radiometer (SMMR)“, also on board Seasat, in- 
dicate that in the future, these instruments may be 
able to estimate surface wind speeds around tropical 
cyclones by observing the sea state (Jones ct aL, 
1979; LipesiT al.. 1979). 

In a previous paper (Kidder al., 1978) we pro- 
posed a microwave technique for estimating surface 
wind speeds at outer radii which has several poten- 
tial advantages over other techniques. First, it is 
suitable for use with the Microwave Sounding Unit 
(MSU) on board the current polar-orbiting opera- 
tional satellites (TIROS-N and NOAA 6). Second, 
it requires only a small amount of computing and no 
special equipment, such as video displays. Third, it 
utilizes routinely available data; no special obser- 
vations are needed. Fourth, the technique is unaf- 
fected by the extensive cloudiness which surrounds 
a tropical cyclone. And finally, it is an objective 
technique based soundly on the radiative transfer, 
hydrostatic and gradient wind equations. Note that 
in related work Grody ct al. (1979) have used micro- 
wave data to estimate weighted mean tropospheric 
wind speeds around Typhoon June (1975). 

In this paper we drive a more exact relationship 
between satellite-observed brightness temperatures 
and surface wind speeds than presented in Kidder 


2 SMMR is also on board the Nimbus 7 satellite. 
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et al. (1978). The theory is first tested using simu- 
lated satellite data constructed from aircraft recon- 
naissance of four hurricanes. The theory is next 
tested using data from the 55.45 GHz channel of the 
Nimbus 6 Scanning Microwave Spectrometer 
(SCAMS) to estimate the radii of 15.4 m s ‘ (30 kt) 
and 25.7 m s ' (50 kt) winds in eight West Pacific 
typhoons during 1975. We call this technique the 
Surface Wind Inference IVom Microwave data 
(SWIM) technique. 

2. Theory 

The high surface winds in a tropical cyclone are 
caused by surface pressure gradients which are in 
turn caused by warm temperature anomalies in the 
upper troposphere, fhe basis of the SWIM tech- 
nique is that microwave radiometers can measure 
the upper level warming which may then be related 
to surface winds. 

(I. Tropical cycloiu' tcmpcralmr .structure 

Shown in Fig. 1 are the mean temperature anoma- 
lies (difference from storm environment) for west 
Pacific typhoons and Ibr Atlantic (West Indies) hur- 
ricanes as composited from rawinsonde data by 
Niinez and Gray (1977). The two main features o*' 
the curves are I) the large positive anomalies peak- 
ing betv/een 25 and 30 kPa and extending several 
hundred kilometers from the storm center, and 2) the 
remarkable similarity in the shape of the profiles 
at different radii. Specifically, the temperature 
anomaly is well approximated by 

7’'(r,z) = «(r)7'(:), (1) 

where 7 (z) is a standard known temperature anomaly 
profile (slightly different in the different oceans) 
and «(/■) is a strength parameter, which will be dif- 
ferent for different storms. That this relationship 
is not an artifact of the compositing procedure may 
be seen in plots of temperature anomalies in in- 
dividual storms (e.g., in LaSeurand Hawkins, 1963: 
Hawkins and Rubsam, 1968; Hawkins and Imbembo, 
1976). 

It must be noted that the profiles shown in Fig. 1 
have been modified slightly from those presented in 
Nunez and Gray (1977). First, becau.se there were 
few statistically significant temperature anomalies 
at low levels, the profiles were extrapolated to zero 
at 100 kPa in the Pacific and at 70 kPa in the Atlan- 
tic. Second, small cold anomalies exist above the 
storms due to overshooting Cb cloud tops and/or 
radiative cooling. Because data are sparse at upper 
levels, this cooling was parameterized as follows. 
In the Atlantic, where storms are shallower 
than in the Pacific, the cooling was assumed to be 
zero at 16.25 kPa and 5 kPa and to be parabolic be- 




Fui. 1. Mean temperature anomalies in V radial bands for 
(a) West Pacitic typhoons and (b) Atlantic (West Indies) hur- 
ricanes (after Nunc/ and Gray. 1^77). 


tween with the amplitude determined by the 20 kPa 
anomaly such that the slope was continuous across 
16.25 kPa. This parameterization fits the data rather 
well (see Kidder, 1979). In the Pacific, the cooling, 
shown schematically with the dotted line in Fig. la, 
was ignored because of statistical uncertainty in 
rawinsonde data at that altitude. It also made no 
difference in the calculations presented below. 
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Fig. 2. Weighting functions for (a) the Nimbus 6 Scanning 
Microwave Spectrometer (SCAMS) and (b) the TIROS-N/NOAA 
Microwave Sounding Unit (MSU) in the 15° North Annual 
Atmosphere. 


b. RcuUoneter characteristics 


The SCAMS instrument is a five-channel radiom- 
eter sensing radiation nominally at 22.235, 31.65, 
52.85, 53.85 and 55.45 GHz (Staelin et al., 1975a). 
The lower two frequencies are used to estimate 
integrated liquid water and water vapor (Grody, 
1976). The upper three channels, which are on the 
wing of an oxygen (O 2 ) absorption band, are used to 
sound the atmosphere. In the absence of scattering 
[due to particles of diameter ^ 0.2 mm (i.e., pre- 
cipitation)] the brightness temperature measured by 
channel /(/ = 3, 4, 5) of the radiometer is given by a 
form of the radiative transfer equation: 


Tn. 


^iTsTi + 


T(.z)Wiiz)dz, 

Jo 


( 2 ) 


where 6/ is the surface emittance, r/ the transmit- 


tance from the surface to the satellite, the surface 
temperature, 7'(z) the atmospheric temperature at 
height z and vi^j(z) a weighting function (Staelin 
c'f al. , 1975a). VK/(z) and e, are both weakly dependent 
on temperature but more strongly dependent on scan 
angle. Fig. 2 shows the weighting functions for the 
SCAMS and for the MSU on board the T1R()S-N 
and NO A A 6 satellites. The 15'' North Annual 
Atmosphere was used to construct these weighting 
functions, and it was assumed that the surface was 
the ocean where e, is 0.5. 

At frequencies which peak in the upper tropo- 
sphere, warm brightness temperature anomalies 
should be observed over tropical cyclones, 'fhese 
warm anomalies were first observed in 55.45 GH/ 
SCAMS data over Typhoon June of 1975 by Rosen- 
kranz and Staelin (1976) and RosenkranzeV al. ( 1978) 
and over other storms by Kidder fV al. (1978). All 
of the above authors concluded that the brightness 
temperature anomalies were not caused by clouds, 
which are far too transparent to cause anomalies of 
several kelvin (Staelin et al.. 1975b; Kidder, 1979) 
or by precipitation, which will be frozen (and thus 
much less able to interact with microwave radia- 
tion) at the levels sensed by the SCAMS 55.45 GHz 
channel. The brightness temperature anomalies are 
directly related to thermodynamic temperature 
anomalies in the storms as shown in Fig. I. The 
object of this paper is to translate ttiese satellite- 
observed brightness temperature anomalies into 
storm parameters of interest to the meteorologist. 

c. Surface pressure ec/uatlon 

Warming in the upper troposphere is responsible 
for the surface pressure drop in tropical cyclones. 
Because satellite-measured brightness temperature 
(at certain frequencies) is a measure of this upper 
warming, brightness temperature ought to be related 
to surface pressure. It was shown in Kidder et al. 
(1978) that 55.45 GHz brightness temperature 
anomalies are statistically correlated with central 
pressure. A more exact relationship may be derived 
as follows. At a frequency for which the transmit- 
tance from the surface to the satellite is vanishingly 
small, the surface term in Eq. (2) may be neglected. 
Thus 


Tn = r nz)W(z)clz. (3) 

Jo 

If the atmospheric temperature is divided into en- 
vironmental temperature T,, plus temperature per- 
turbation T , and if Eq. (1) is used to express T', 
we can write 


Tn 


T,(z) W(z)dz + 
Jo 


" T'{r,z)W(z)dz 
Jo 


(4) 
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or 

7 „ - 7 „, t «{r) I i\z)W(z)dz, ( 5 ) 
Jo 

where is the environmental brightness tempera- 
ture. In the above, the small variation of the weight- 
ing function with temperature has been neglected. 

The surface pressure outside of the radius of maxi- 
mum wind can be obtained from the hydrostatic 
equation.’^ If the hydrostatic equation is integrated 
from the surface to a height such as the 5 kPa level 
which is undisturbed by the storm (Frank, 1977), 
we have 

\nps - In/^/ + — . (6) 

R Jo T(z) 


wSubstituting for 7'(z) as above and using the binomial 
expansion, keeping only terms linear in «(;*), yields 


\nps 


In/;,. - a(r) ^ 


f(z) , 

clz 

,0 7,.“(z) 


(7) 


where is the environmental surface pressure, and 
Hr the height of the 5 kPa surface. The unknown 
tx{r) may now be eliminated between Eqs. (5) and 
(7) to yield 

A In/;, ^ (8) 

where A indicates a departure from the environ- 
mental value and A is given by 

iLrjQiL,, 


f(z)W(.z)dz 

. 0 

F,q. (8) relates surface pressure anomalies to 
anomalies in a satellite-observed quantity (bright- 
ness temperature). Its usefulness depends on the 
characteristics of the coefficients. Of chief concern 
is the uncertainty in the calculated value of A. An 
alternate method for calculating A which is slightly 
more accurate than Eq. (9) [because it takes into 
account the small variation of W(z) with tempera- 
ture] is to calculate the brightness temperature for 
a storm environment, then to perturb the environ- 
mental temperature with one of the profiles in Fig. 1 
and calculate the surface pressure and brightness 
temperature anomalies. When this procedure was 
carried out using Pacific and Atlantic environmental 
soundings from the work of Nunez and Gray (1977) 
and atmospheric absorption coefficients calculated 
using the formulas of Liebe and Gimmestad (1978) 
and Liebe ct tiL (1977), the values of A listed in 


vSee Kidder (1979) for a demonstration that on space scales of 
100 km the hydrostatic equation is valid in tropical storms. 


Taiui I. Values of the coefficient /t (K ‘) as a function 


of radius. 


Radial band 


(deg latitude) 

55.45 (HI/ 

54.96 (HI/ 

Atlantic 




0^1 

1.13 X 

10 “ 

0 66 X 10 - 


1.10 


0.66 

2^3 

1.07 


0.64 

.3-4 

1.03 


0.62 

4^5 

1. 01 


0.62 

5 6 

0.94 


0..M 

6-7 

0.85 


0.51 

Mean 

1.02 X 

10 - 

0.61 X 10 - 

Standard deviation 

0.09 X 

10 “ 

0.06 X 10 - 

Pacific 




0^1 

0.90 X 

10 “ 

0.82 X 10 - 

!^2 

0.94 


0.83 

2-3 

0.98 


0.84 

3-4 

0.99 


0.86 

4^5 

1.02 


0.87 

5^6 

0.95 


0.85 


0.83 


0.80 

Mean 

0,95 X 

to “ 

0.84 X 10 - 

Standard deviation 

0.06 X 

10 " 

0.02 X 10 “ 


Table 1 resulted. There seems to be no systematic 
variation in A with radius, and the standard devia- 
tions of the samples are less than I09f of the mean 
values. Furthermore. A seems to be fairly insensi- 
tive to changes in the environment. Changing to the 
U.S. Standard Atmosphere changes A only ^ 
Similarly, switching an Atlantic temperature anomaly 
profile to a Pacific profile only changes A by 11%. 
It thus appears that the coefficient A can be cal- 
culated from existing data, and it has a rather small 
uncertainty. 

Kidder et al, (1978) showed a high correlation 
(-0.86) between central surface pressure and the 
brightness temperature anomaly at the center of the 
storm. Although there are problems with this rela- 
tionship due to the eye having a variable diameter, 
the high correlation is explained by Eq. (8). 

c/. Surface wind speed 

If one has an estimate of the surface pressure 
gradient, one ought to be able to estimate the sur- 
face wind speed. Our simple approach to this prob- 
lem is the following. Although surface pressure in 
tropical cyclones is neariy symmetric about the cen- 
ter(Frank, 1977), wind speeds are notoriously asym- 
metric. Most of that asymmetry at radii V\ how- 
ever, can be explained by the motion of the storm 
(George and Gray, 1976). Therefore, we assumed 
that at gradient level ( 85 kPa) the relative winds 
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(Storm motion subtracted) were symmetric and to a 
close degree in gradient balance,' i.e,. 



r 


lUr, 


() In/)., 


( 10 ) 


where V',, is the relative wind speed at gradient level. 
r the radius, R the gas constant. 7,; the temperature 
at gradient level,/the Coriolis parameter, and /),, the 
surface pressure. The surface pressure has been sub- 
stituted for the pressure at gradient level because, 
by the hydrostatic equation, the radial gradient of 
the logarithms of the two pressures will be identical 
if the radial gradient of temperature between the two 
levels is small, which is the case (I'ig. 1). .Substitut- 
ing from Hq. (8) yields 



r 


+ fV, 


ARTu 


Or 


( 11 ) 


Kq. (11) relates relative wind speed at gradient level to 
satellite-observed brightness temperature gradients. 

If satellite data were noise-free, one could take 
any two brightness temperature observations near 
a tropical cyclone and estimate the gradient level 
wind speed between them from Hq. (11). Unfor- 
tunately, the SCAMS instrument has about tO.5 K 
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noise equivalent brightness temperature which would 
result in e.xtremely noisy winds. It is therefore 
desirable to smooth the data before calculating the 
winds. One way to accomplish this smoothing is to 
assume a functional form for wind speed with radius. 
A suitably simple function proposed by Hughes 
(I9.S2), Riehl (ly.M, 19W). and others is 

VO'. (12) 

The exponent .v would be 1 for parcels conserving 
their angular momentum. Near the surface where 
friction is large, v is nearO..*' (Riehl. I%.V Oray and 
.Shea, l‘)7.1). The effect of changing ,v will be dis- 
cussed later. InscrHng Hq. (12) in Hq. (11) and in- 
tegrating with respect to r holding ( \ A and /,, 
constant, gives 

/« lARi,.) '(c-~- f /,. (1.1) 

' 2.V 1 -V 

where is an integration constant, (liven a lieki 
of satellite-measured brightness temperature around 
a tropical cyclone, one can average the temperatures 
in radial bands about the storm center (an additional 
smoothing technique) and calculate the which 
gives the best least squares lit to the data. This (' 
is the real, positive rook’ of 

_L (V /• 

N 



= . -(V.VHV 
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0. (14) 


where y =ARTaTi,. 

One additional parameter is necessary. We used a 
simple shearing parameter to relate the surface wind 
speed to the wind speed at gradient level, i.e., 

V., = /xV„. (l-S) 

Based on the work of Bates (1977) and Gray and 
Frank (1978), 0.7 was chosen as the value of fx. Be- 
cause wind speeds above gradient level change slowly 
in comparison with those below gradient level 
(Frank, 1977; Bates, 1977), accurate knowledge of 
the height of the gradient level is not required. 

The SWIM technique for estimating surface wind 
speeds from satellite-observed brightness tempera- 
tures can then be summarized as follows: 

1) Average the brightness temperatures in radial 
bands about the storm center. [We used 56 km ('/i“) 
bands.] 

2) Choose an appropriate value of .v (see next 
section). 

Tlii.s .issumption i.s not true at the surface where friction is 
large. .Surface wind speeds will be derived from gradient level 
wind speeds by a simple parameterization of frictional effects. 


3) Calculate C with Hq. (14). 

4) Calculate surface wind speeds using Hqs. (12) 
and (15). 

5) Add the .storm motion to the relative winds to 
obtain the asymmetric wind field. 

The technique was tested using, first, simulated 
satellite data constructed from temperature, pres- 
sure and height data measured by hurricane recon- 
naissance aircraft, and second, using 55.45 GHz 
SCAMS data. 

3. Results 
a. Aircraft data 

Although satellite data are prevalent over tropi- 
cal cyclones, one invariably finds that the number 
of cases in which one has both satellite data and 
ground truth data is extremely small. Also ground 


It is observed that Hq. (14) has only one real, positive root, 
the other two are generally real and negative, although they arc 
sometimes complex conjugates. Occasiontilly for weak storms, 
no positive real root is found. 
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Tabu 2. Summary of aircraft reconnaissance data 
from hurricanes used in this study. 




Flight levels 

Intensity 

Storm 

Date 

(kPa) 

change 

Cleo 

18 Aug 1958 

80. .S6. 24 

Steady 

Helene 

26 Sep 1958 

80. .S6, 25 

Deepening 

Hilda 

1 Oct 1964 

90. 75, 65. 50. 18 

Deepening 

Hilda 

2 Oct 1964 

90. 70. 65. 20 

Filling 

Inez 

28 Sep 1966 

95. 75. 65. 50, 20 

Deepening 


resolution of microwave data is rather coarse — 
145 km at nadir for SCAMS, for example. F’or these 
reasons, we decided to first test the SWIM technique 
using simulated satellite data which had good ground 
truth data and high spatial resolution. Aircraft recon- 
naissance data collected for the National Hurricane 
Research Project (Gray and Shea, 1976) were utili/ed 
for this simulation. The data were screened for the 
following criteria: 1) There had to be at least three 
flight levels; 2) one of the flights had to be near the 
level of maximum warming (25-30 kPa): and 3) one 
of the flight levels had to be near the surface to 
measure low-level winds. Between 1957 and 1969 
there were five days (four storms) on which such 
data were taken (Table 2). 

The satellite data were simulated as follows. At 
each 4.6 km (2.5 n mi) radial grid point, smoothed 



Fig. 3. Calculated 55.45 GHz brightness temperatures (circles) 
and best fit line for four hurricanes. 



10 20 30 40 60 80 100 

RADIUS (km) 

Fig. 4. Calculated gradient level ( 85 kPa) wind speeds (line) 
and observed low-level wind speeds (circles) for four hurricanes 
(note that both scales are logarithmic). 

vortex-averaged values of temperature, pressure 
and height (at the lowest flight level) from Gray and 
Shea (1976) and relative humidities from the mean 
hurricane environment (Kidder, 1979) were used to 
construct soundings. Between flight levels, the 
equivalent potential temperature was assumed to 
vary linearly with pressure (LaSeur and Hawkins, 
1963), and above the uppermost flight level, the 
mean environmental sounding was used. Finally, 
55,45 GHz brightness temperatures were calculated, 
from Eq. (2), for each grid point. The result was 
simulated, vortex-avei aged, 55.45 GHz satellite data 
with a spatial resolution of 4.6 km (2.5 n mi). 

Simulated data outside of the radii of maximum 
wind were inserted into Eq. (14) to calculate C for 
a: = 0.5 and A = 1.02 X lO"- K“^ The offset T, 
was also calculated, and the simulated brightness 
temperatures along with the best-fit line from Eq. 
(13) are shown in Fig. 3. It can be seen, particularly 
in the case of Inez, that Eq. (13) fits the data well. 
The wind speeds at gradient level (-85 kPa) were 
then calculated with Eq. (12), and they are com- 
pared with the observed, vortex-averaged wind 
speed at the lowest flight level (Fig. 4). With the 
exception of Hilda on 2 October 1964, the agreement 
between observed and calculated wind speeds is 
excellent with rms errors of 2-3 m s ^ 
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Tabi i 3. Typhoons used in this study (all 1975). 


Central 
location vU 

Date/time Minimum time of peak 


Peak (GMT) sea level intensity 

intensity of peak pressure , 


Name 

(m s ') 

intensity 

(kPa) 

(’N) 

V'li) 

June 

82 

19 Nov/ 1 200 

87.5 

13.2 

14 1. 0 

Phyllis 

62 

14 Aug/ 1800 

92.5 

24.1 

1.37. 1 

Rita 

41 

22 Aug/1200 


32.9 

134.4 

Tess 

49 

4 .Sep/1800 

94.5 

23.0 

147.6 

Winnie 

.33 

10 Sep/0600 


31.0 

162.8 

Alice 

39 

17 Sep/1200 

97.3 

15.4 

1.32.1 

Betty 

49 

22 Sep/0000 

94.7 

22,6 

123.6 

Cora 

54 

4 Oct/1800 

94.3 

30.3 

133.2 


The case of Hilda on 2 October 1%4 is interesting 
because it is the only filling storm. It is possible that 
storm temperature structure is different enough in 
later stages that A should have a different value. 
It is also possible that the process of constructing 
soundings mentioned above failed to detect the true 
magnitude of the upper level temperature anomaly. 
This may have caused the calculated brightness 
temperature gradient, and thus the calculated winds, 
to be too low. More research is needed to adequately 
resolve this problem. 

In summary, it appears that with 4-5 km resolu- 
tion, low-noise satellite data, one could accurately 
deduce low-level tropical cyclone wind speeds. We 
next attempted to estimate wind speeds using Nim- 
bus 6 SCAMS data. 

/?. Nimbus 6 data 

In ourprevious study (Kidder 1978), we had 
collected 1 16 Nimbus 6 SCAMS images of hurricanes 
and typhoons during 1975. These images were 
screened to meet two criteria: 1) the scan angle to 
the storm center had to be less than or equal to 
21.6'', and 2) an independent estimate of the radius of 
15.4 m s ‘ (30 kt) wind had to be available. Only 20 
images met both criteria, and only one ofthe.se was 
in the Atlantic; the remainder were Pacific typhoons. 
To make the data set more uniform. Tropical Storm 
Gladys (the only Atlantic storm) was eliminated, 
leaving 19 images of eight typhoons (Table 3). This 
is a rather small data set but it is large enough to 
indicate the potential of the SWIM technique. 

Before calculation of the wind speeds, the 55.45 
GHz brightness temperatures were corrected for 
scan angle with an additive correction developed by 
Rosenkranz et al. (1978), The center of the storm 
was cho.scn as the maximum brightness temperature 
within 145 km (one scan spot) of the interpolated 
best track center (1975 Annual Typhoon Report). 
The brightness temperatures were then aziniuthially 
averaged in 56 km (0.5'' latitude) bands from 1 1 1 km 
( 1°) to 778 km (7°). The brightness temperature in the 


very center of the storm was not used in the fitting 
process. lup ( 14) was used to calculate C’ for a 0.5. 
A 0.95 X 10 ‘ K * and 17.2T. Finally, the 
radiiof 15.4 ms * (30 kt) and 25.7 ms M.5()kt) winds 
were calculated using /x 0.7 to estimate surface 
wind speeds. 'Fhese calculations are compareti with 
observations in Fig. 5. 

The observed radii were taken from logs kept at 
the Joint Typhoon Warning Center on Guam. Hach 
6 h the forecaster on duty estimates the radii of 30, 
50 and 100 kt winds using all available data. There 
arc two problems with these observations. First, 
they are not coincident in time with the satellite 
image. The observational values plotted in big. 5 
have been averaged both around the storm and over 
the two 6 h estimates between which the satellite 
observation was made. Second, the observations are 
subjective estimates which introduces some variance 
and a bias. Because safety is the prime concern, the 
observations tend to be overestimates, fhe average 
value of this error is unknown. 




5. Observed versus estimated average radius ol'(ii) 15.4 m 
s ‘ (30 kt) and (b) 25.7 ms* (50 kt) winds in eight typhoons dur- 
ing 1975 for A 0.5. 
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The regression lines in Fig. 5 explain 63 and 749? 
of the variance for 25.7 m s ' (50 kt) and 15.4 m s ' 
(30 kt) winds, respectively. The standard errors of 
estimate are 68 and 93 km, respectively, which are 
reasonably good, and which may be influenced by 
the nature of the observed values. The intercepts 
are rather large, but some of this error is accounted 
for by the natural overestimation in the observed 
radii discussed above. 

If the theory is correct, the slopes of the regres- 
sion lines should be near 1. The fact that the slope 
for the radius of 15.4 m s ' wind is somewhat less 
than 1 leads to speculation on the best value for the 
exponent .V in the wind speed profile. Within 100 km 
of the storm center, 0.5 is a reasonable value for.v, 
but Malkus and Riehl (1959) suggest that between 
200 and 500 km .v should be 0.6 to 0.7. Even though 
very few data points were available, an analysis was 
done to determine the effects of changing .v. 

Table 4 shows the square of the regression coef- 
ficient, the standard error of estimate, the slope and 
the intercept of the regression lines as functions of.v. 
In general, increasing .v slightly increases the frac- 
tion of explained variance, increases the slope, de- 
creases the intercept and slightly decreases the 
standard error of estimate. For the radius of 25.7 m 
s ' winds, an v in the range 0.4-0.5 produces the 
desired slope near 1.0. For the radius of 15.4 m s' ' 
winds, however, it appears that an.v in the range 0.6 
to 0.7 is a better choice. Fig. 6 shows the effect of 
changing .v to 0.7. A better fit seems to be achieved 
for the radius of 15.4 m s ' winds, but the opposite 
is true for the radius of 25.7 m s ' winds. It appears, 
then, that for calculating the radius of 15.4 m s ' 
(30 kt) winds by the SWIM technique one should 
u.se .V = 0.7, but for the radius of 25.7 ms' (50 kt) 
winds one should use .v = 0.5. 

The 109? uncertainty in A causes approximately 
a 6-7% uncertainty in the wind speeds. The result- 
ant uncertainty in the radius of a particular wind 


Table 4. Variation of regression statistics as a function of.v. 


,v 


Intercept 

(km) 

& 

(km) 

Slope 


Radius of 15,4 m ,v" ‘ 

[50 kt) wind 


0.4 

0.680 

349 

103 

0.551 

0.5 

0.736 

317 

93 

0,737 

0.6 

0.771 

285 

86 

0.915 

0.7 

0.801 

251 

81 

1.090 

0.8 

0.820 

217 

77 

1.265 


Radius of 25,7 m a ’ 

(50 kt) wind 


0.4 

0.579 

111 

72 

1.058 

0.5 

0.631 

92 

68 

1.107 

0.6 

0.663 

74 

65 

1.165 

0.7 

0.685 

55 

63 

1.231 

0.8 

0.701 

36 

61 

1.305 




Fici, 6, As in Fig. .S except for.v - 0.7. 


speed is the latter figure divided by .v, i.e., a 139? un- 
certainty for .V = 0.5 and a 10% uncertainty for .v 
= 0.7. A better knowledge of the coefficient A would 
allow more accurate determinations of wind speed. 

4. Conclusions 

From Figs. 5 and 6, it appears that the SWIM 
technique is currently capable of estimating the radius 
of 15.4 m s“' (30 kt) wind to within about ±80 km and 
the radius of 25.7 m s“' (50 kt) wind to within about 
±70 km. These results are promising, and indicate 
that further research is required. 

The SWIM technique should be tried using 
TIROS-N and NOAA 6 data for several reasons. 
First, the MSU is a more sensitive instrument for 
estimating surface pressure grijdients than was the 
.'^CAMS. The peak of the 54.96 GHz channel weight- 
ing function is closer to the peak of the temperature 
anomaly than the peak of the 55.45 GHz SCAMS 
weighting function (Fig. 2). Also, the 57.95 GHz 
channel could be used to measure the cooling above 
the storm. A two-parameter estimation of the sur- 
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face pressure anomaly of the form 

A \r,p, AA7„, + BM„, (16) 

should be more accurate than a single-parameter 
estimation. Second, since the launch of NOAA 6 
there are four images per day of tropical cyclones, 
about half of which pass closely enough to the storm 
center on the average to apply the SWIM technique. 
And finally, the TIROS-N.'NOAA series are the cur- 
rent operational satellites from which microwave 
data are routinely processed in near real time. This 
technique is also suitable for use with the microwave 
sounder on hoaid the new DMSP satellite, and further 
te.sting may be done with those data. 

In summary, we think that the microwave tech- 
nique described in this paper, which we have called 
the SWIM technique, has the potential to become a 
valuable tool at the di.sposal of the forecaster. 
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ABSTRACT 

An analysis of tropical cyclone forecast track errors shows that the largest errors are typically 
associated with storms undergoing turning motion. This paper presents results obtained from a com- 
posite study of tropical cyclones occurring in the West Indies during I%I-77. Storms which under- 
went a left or right turn or moved straight for a period of at least 36 h were studied. Wind ields at 
different levels in the atmosphere around these storms were investigated 

When a storm begins to turn, the environmental How at 500 mb and the average tropospheric wind 
(between 200 and 900 mb) around the cyclone at 5-ll‘-* latitude radius is cyclonic for a left turning 
and anticyclonic for a right turning storm. At 24-36 h before a storm makes a left turn, there exists a 
large positive vertical wind shear around the cyclone between the upper (200 mb) and lower (900 mb) 
troposphere in the direction of storm motion, while the opposite occurs with a right-turning storm. 
Straight-moving cyclones do not show such a shear pattern. Statistical tests show that these results are 
signilicant at the 950' level or higher. 

Tropospheric mean temperature fields around 13 tropical cyclone turning cases in the Atlantic and 
Pacific Oceans derived from the Nimbus 6 microwave sounder data during 1975 were also studied. 
Temperature gradients across these storms indicate (through the thermal wind relationship) vertical 
wind shear profiles similar to those obtained from the composite. 

These results suggest or verify previous ideas that 1) by measuring certain parameters around a storm 
(sense of surrounding wind rotation, vertical wind shear between 200 and 900 mb, or gradient of 
tropospheric mean temperature) one may be able to make a better 24-36 h forecast of cyclone 
turning motion; 2) the turning motion of tropical cyclones is controlled by large-scale flow fields 
surrounding them; and 3) there seems to be a time lag between the changes in the environment and 
the response of the storm center to such changes. 


1. Introduction 

To predict the :ig motion of tropical cyclones 
24-36 h in advance is very difficult and quite often 
unsuccessful. Statistical prediction schemes (e.g., 
Hope and Neumann, 1970; Neumann and Lawrence, 
1975) have great problems in forecasting such turn- 
ing motion because of their strong emphasis on 
climatology and persistence. Numerical models also 
fail to predict cyclone turning motion most of the 
time. As a result, the largest tropical cyclone track 
forecast errors* are usually associated with storms 
undergoing a turn. 

A detailed study of the characteristics of the sur- 
rounding environment of cyclones undergoing turn- 
ing motion may give insight into ways of devising 
new schemes which may reduce forecast errors 
associated with turning motion. This paper presents 

' Present affiliation: Laboratory for Atmospheric Research, 
University of Illinois at Urbana-Champaign, Urbana 61801. 

' Forecast error is defined as the deviation of the forecast 
position of a tropical cyclone from its best-track verifying 
position. 


information on the environmental flow patterns 
around tropical cyclones prior to and at the time 
when these storms undergo turning motion. It will 
be shown that significant differences in the large- 
scale surrounding wind fields exist at 500 mb, and 
also at 200 and 900 mb levels between left-turning, 
straight-moving and right-turning storms. Tempera- 
ture sounding data from the Nimbus 6 Scanning 
Microwave Spectrometer also provide evidence of 
some of these differences. 

2, Forecast errors associated with turning motion 

To examine this forecast problem more closely, 
a special study of tropical cyclone forecast errors 
occurring in the West Indies during the period 
1961-77 was made. 


a. Selection of cyclone tracks 

The data sample consists of tropical cyclones 
in the West Indies west of 55°W with maxi- 
mum sustained winds ^ 18 m s~*. The direction of 
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Fig. 1. Tracks of tropical cyclones used in this study which made a left turn. The asterisk on each track indicates the turn time 
T. The solid circle is the 0000 GMT position (with the date next to it) and the open circle is the 1200 GMT position. 


motion (D) of each storm at a standard time T in 
hours (0000 or 1200 GMT) is defined as D(T). Three 
s turn classes were then chosen according to the fol- 

lowing criteria: 

• Left-turning: ,D(T + 12) - D(T) < -20° 

1 Straight-moving: -10° < D(T + 12) ~ D(T - 12) 

, < 10 ° 

^ Right-turning: D(T + 12) - D(T) > 20°. 

f The time T was then defined as the turn time. In 

this way, 16 left-turning, 33 straight-moving and 28 
right-turning cases were chosen. Tracks of these 
storms are shown in Figs. 1, 2 and 3. 

* /;. Forecast errors of sample storms 
s 

1 The 24 h National Hurricane Center (NHC) 

official forecast errors issued for each storm at three 
time periods [24 h before turn time (T “■ 24), the 
turn time (T) and 24 h after turn time (T + 24)] were 
analyzed. Table 1 gives the average result for each 
> turn class. The ‘'special right turn" category con- 

sisted of cases where the average 24 h forecast 
f error was greater than 350 km. It can be seen that a 


considerable jump in the forecast error occurs at 
turn time in comparison with the forecast error for 
straight-moving cyclones. The forecast errors are 
even larger for 16 of the 22 special right-turn j 

cases, jumping from 148 km before the turn to 417 | 

km at turn time. 

Fig. 4 shows the scatter of 24 h foreca.>t posi- 
tions made at turn time and the corresponding \ 

verifying positions for each turn class. Table 2 | 

Table I. Average 24 h oflicial tropical cyclone track forecast 
errors (km) issued by the National Hurricane Center, Miami. T 
Is the time when the storm starts to turn. See text for further 
explanation of turn clas.sification, Note that the number of 
analyzed cases are different from those mentioned in Section 2a 
because forecast errors for some cases are not available. 


24-h Forecast error (km) 
for forecasts issued at 
the following times 


Turn classification 

T - 24 

T 

T + 24 

Left turn ( 10 cases) 

235 

289 

206 

Straight (23 cases) 

148 


196 

Right turn (22 cases) 

178 

324 

239 

Speciiil right turn ( 16 cases) 

148 

417 

245 
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Fia. 2. As in Fig. I except for straight-moving tropical cyclones. 


gives the average directional deviation of these 
forecast positions from the verifying positions and 
the extrapolated tracks. The forecast positions for 
right-turning storms are nearly always to the left of 
the actual storm track [Directional Deviation (DD) 

- -38'"]. They also tend to cluster around the 
exti’apolated position of the storm from its cun-ent 
path(DD = 2^^). The strong bias based on persistence 
is evident. All forecast positions for left-tuiming 
storms arc to the right of the actual ti’ack (DD 

- 50*^). It is interesting to note that, on the average, 
the forecast actually calls for a right-turn (DD from 

Taulu 2. Directional deviation (DD) of the r. .m 24 h forecast 
position made at turn time from the mean verifying position and 
the mean extrapolated track. A positive number means the fore- 
cast position is to the right of the verifying position or the 
extrapolated track. 


Turn 

classilicatjon 

Directional deviation of foreca*«t (deg) 

l^'rom mean 
extrapolated (rack 

From mean 
verifying position 

Left turn 

21 

50 

Straight 

9 

11 

Right turn 

2 



extrapolated track is 2T). For straight-moving 
storms the forecast, verifying and extrapolated posi- 
tions are nearly the same. Although there are natural 
limitations on the predictability of storm movement, 
it seems obv ious that the skill of forecasting turning 
motion is not satisfactory. This is probably a result 
of insuflicient observational knowledge of the char- 
acteristics of the cyclone's surrounding environ- 
mental How patterns and their changes. It is apparent 
that more research on cyclone turning motion is 
needed. 

3. Data compositing 

Because of the scarcity of oceanic rawinsonde 
data, representative information on the characteris 
tics of surrounding flow patterns of individual tropi- 
cal systems is usually very deficient. To overcome 
this severe data problem, a rawinsonde compositing 
approach was employed. This averages data from 
storm systems having similar turning characteristics. 
The rationale behin(J this type of data compositing 
has been discussed in a number of recent papers (e.g., 
Williams and Gray, 1973; Frank, 1977; and others). 
Composited data samples enable reliable quantita- 
tive analyses to be made, and allow meaningful 
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Fig, 3. As in Fig. I except for right-turning tropkal cyclones. 


conclusions to be drawn about the '‘average" differ* tween turning classes (if present) will give hints to 
ences in the surrounding environment between dif- forecasters as to which surrounding cyclone param- 
ferent data sets. Compositing tends to suppress eters to monitor in individual cases, 
random data noise and isolates mean characteristics For each turn classification in wSection 2a rawin* 
common to each composite class. Differences be- sonde data at the following time periods wer; 


LEFT TURN 


STRAIGHT MOVING 


RK3HT TURN 


T^tj-FORECAST 
• F^?FSENT 

T+ 24L*^ OBSERVED 




”°o *^^+ 24 (g^ ''♦^OBSERVED TRAcT 


T+24 4 

4 ?K^0eSERVED 

°.V- :>>,,Pf>6SENT 

4 4**/ T TftACK 

. •forecast * 


5' 4' 3* 2* i* O’ 

DISTANCE IN •UTITUDE 


10’ 9’ 8’ 7’ 6’ 5’ 4* 3* 2’ J* O’ 

DISTANCE IN “LATITUDE 


5* 4’ 3* 2* I’ 0* 
DISTANCE IN “LATITUDE 


Fig. 4. National Hurricane Center 24 h forecast positions (denoted by +) made at the turn lime for individual storms in 
each class and their actual positions (denoted by open circles) at the lime T ^ 24. The is the mean forecast position. The 
dashed arrow is the 24 h mean forecast track. The solid arrow from T to T ^ 24 is the mean distance (and direction) the 
storm has traveled during the 24 h. The dotted arrow indicates the average displacement of the forecast position from the 
actual position. 
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composited; when the storm starts to turn (turn 
time, T), and 12, 24 and 36 h before turn time 
(T - 12, T - 24 and T - 36 respectively), as shown 
schematically in Fig. 5. R:i\/insondc stations used in 
the data compositing are shown in Fig. 5. 

Compositing is performed on a 15"’ latitude radius- 
cylindrical grid consisting of eight octants and eight 
radial bands, as shown m Fig. 7. The center of the 
grid corresponds to the storm center. The grid is 
rotated so that octant 1 is aligned along the direc- 
tion of storm motion. See George and Gray (1976) 
for a more detailed description. 


Fig. 5. Idealized picture of the three turn classes of tropical 

cyclones and of the time periods prior to turn which were = Hereafter, all distances are referred to in degrees Latitude 
composited. (I°latitude= 111. 1 km). 



Fig, 6. Northwest Atlantic rawinsonde data network. 
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Fici. 7. Grid used for compositing rawinsonde data. The 
arrow points in the direction of storm motion. Outer numbers 
denote octants. Numbers inside grid indicate distances from the 
center in degrees latitude. 


Each wind vector is resolved into two components, 
one parallel to the direction of storm motion (V/.) 
and one normal to it (Vv)i tts shown in Fig. 8. All 
wind observations falling within a particular octant 
and radial band are averaged to give the composite 
and Vs within that sector. These types of 
rawinsonde averaging processes are described in 
greater detail by Frank (1977). 

4. Surrounding How at turn time 

It has been known for many years that tropical 
cyclone motion is well related to mid-tropospheric 
surrounding wind patterns. In a composite study of 
tropical cyclone motion stratified according to direc- 
tion, latitude, speed, intensity and intensity change, 
George and Gray (1976) found that 500 mb is the best 
steering level for direction and 700 mb is the best foi 


Vp 



STORM MOTION SYSTEM 

MOTION COMPONENT VECTORS 

VECTOR 

Fui. 8. Parallel (\'/.) and perpendicular (V\) components of 
a wind vector showing their relation to th^ storm motion 
vector. 

cyclone speed. To better refine the present knowl- 
edge of the steering flow associated with cyclones 
undergoing turning motion, composite flow fields 
at 900, 700, 500 and 200 mb at turn time were 
analyzed. Although 500 mb is found to be the best 
steering level, flow fields at the upper (200 mb) and 
lower (900 mb) tropospheric levels also appear quite 
useful. 

a. Flow field at 500 mh 

Fig. 9 shows the composited 500 mb streamline 
patterns at turn time for the three turning classes. 
The flow field ahead of the storm can be seen to be a 
very good indicator of current directional change. 

The Vs components of the winds at 500 mb in the 
front and back octants are shown in Fig. 10. The 
number next to each arrow indicates the magnitude 
of the Vs component (m s"M averaged between 
5-ir radius. For each classification, if a weighted 
average of V\ in the front three and back three octants 
[Vs in octants 1 (front) and 5 (back) are given a 



Fig. 9. 500 mb streamlines for the three turn classes at turn lime. Open arrows indicate the instantaneous direction of storm motion. 
Solid arrows indicate the movement of the storm during the next 12 h. 


! 

I 
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(FRONT AVERAGE) 
-(BACK AVERAGE) “15 



Fici* 10. Sum of Vs (wind component (m s ‘) normal to direction of storm motion] at 500 mb within 
5-1 1" radii from the storm center in various octants at the turn time T. The length of the arrow is 
proportional to the magnitude of the wind. See text for description of the averaging process. 
Open arrows indicate the instantaneous direction of storm motion and solid arrows the movement 
of the storm during the next 12 h. 


weight of 1 and those in the right and left octants a 
weight of 0.5] is calculated, the front octant aver- 
age minus the back octant average will give a hori- 
zontal V,v shear value. This is shown at the bottom 
of the figure. These shear values are significantly 
different between the three turn classes and they 
well infer turning motion at turn time, 

With decreasing oceanic middle level ravvinsonde 
and propeller aircraft data, such 500 mb level 
information on individual cyclone cases is usually 
not available. Since upper and lower tropospheric 
wind data are becoming more available with time 
from satellite and jet aircraft observations it may be 
more useful to try to establish steering fiow relation- 
ships at upper and lower tropospheric levels. 

h. Upper and lower tropospheric flow 

The 5- ir, radius average of the fronc and back 
octants V.v (200 mb) and Vx (900 mb) for the three 
turn classes at turn time is shown in Fig. II. The 
differences in horizontal Vx shear values between 


turn classes follow the same trend as those at 
500 mb. Thus, it seems possible to use the average 
200 and 900 mb winds to infer turning motion at 
turn time. 

5. Surrounding flow before turn time 

A study of the composite Vx flow fields 24 or 
36 h before turn time at various levels reveals no 
significant differences between the turn classes. The 
Vp component of the wind at 500 mb and the upper 
(200 mb) and lower (900 mb) troposphere, however, 
shows important relationships with cyclone turning 
motion 24-36 h before it takes place. 

a. 500 nth flow 

The Vp components of the winds (parallel to the 
direction of storm motion) at 500 mb in the front 
left (2), back left (4), back right (6) and front right 
(8) octants are found to best correlate with turning 
motion 24-36 h before turn time. These components 


LEFT STRAIGHT RIGHT 





(FRONT AVERAGE) ^ 

-(BACK AVERAGE) “5 +13 


+22 


Fig. II. As in Fig. 10 except V\ is the sum of the average winds between 200 
and 900 mb at the turn time (T). 


JUNi: 1980 


CHAN. GRAY AND KIDDKR 


785 


LEFT STRAIGHT RIGHT 






Fig, 12. Sum of 500 mb Vf. [wind component (m s ') parallel to direction of storm motion] 
within 5- 1 T radii from the storm center in octants 2, 4, 6 and 8 for the three turn classes at 24 h 
(T - 24) and 36 h (T - 36) before turn time. The heavy arrow indicates the direction of storm 
motion. See text for description of calculating the front and back shear. 


are shown in Fig. 12. At 24 h before turn time (T 

- 24), there exists at the front of the storm a strong 
cyclonic shear of Vp [F/, (right) - F/)(left)] for the 
left-turning storms. This shear is denoted by Wp! 
A//, where A// represents the distance between 
octants 2 and 8 or 4 and 6. This front Vp shear is 
about zero for the other two turn classes. Behind 
the cyclone, the right minus left Vp shear is cyclonic 
for all three classes. The magnitude of this Vp shear 
is largest for right-turning storms. Subtracting the 
Vp right minus left shear behind the storm from that 
in front of the storm, that is, [{V p^ - Vp.^i\n 

- (Vp,, - K/.,)/A//], where the numeric subscripts 
denote octant numbers, a front minus back Vp shear 
can be obtained which might be used to identify 
turning motion 24 h before it takes place: a positive 
value for left-turning, small negative for straight- 
moving and large negative value for right-turning. 

At 36 h before turn time (T - 36) a similar differ- 
ence in the front minus back Vp shear exists between 
the three turn classes, as indicated in the lower half 
of Fig. 12. Notice that the front minus back V/* shear 
for straight-moving storms at T 36 is the same as 
that for left-turning storms at T - 24. This appears 
to reduce the distinguishing power of this shear to 
forecast turning. However, given the current low 
skill of forecasting turning motion, this 500 mb shear 
information may still be useful. 

If aircraft reconnaissance flights were available 
to better measure these middle-level outer radius 


shear patterns, then an improvement in forecasting 
turning motion might result. Such aircraft recon- 
naissance requires special propeller aircraft capa- 
bility and will likely be unavailable in many or most 
cyclone forecast situations. By contrast, satellite- 
derived winds in the upper and lower troposphere 
should be more easily available. 

h. Upper and hnver tropospheric flow 

The component of the 200-900 mb vertical wind 
shear parallel to the direction of storm motion 
24-36 h before turn time is significantly different 
between the three turn categories. 

Fig. 13 shows the average Vp wind within 7-11° 
radius** from the storm center in octants 1 (front) 
and 5 (back) at 200 and 900 mb 24 h before the turn 
time (T ~ 24). In front of the storm, right-turning 
cyclones show a large negative wind shear between 
these two levels [F/>(2()0 mb)-F/>(900 mb)] while 
both left-turning and straight-moving storms have 
positive wind shear values. Behind the storm, the 
vertical wind shear is negative for straight-moving 


^ The reasons for using 7~ 1 1“ are twofokl: I) since this fiow is 
related to directional change of the storm 24-36 h ahead, one 
would expect this relation to he more significant further out from 
the storm cen'er, where the environmental steering flow in- 
fluences are much more dominant than that of the storm, and 
2) the value ofVVin this radial band correlates best with the direc- 
tional change of the storm. 
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Fici. 13. Cross-sectioniil view of the average V/. wind components im s ') at 200 and 900 mb 
within 7-U'^ radii from the storm center in octants 1 and 5 24 h before turn time iT - 24). Ihe 
circled numbers are the values of the vertical wind shear lV'/.(200 mb)- '0 mb)]. 


storms and slightly positive for the other two turn 
classes. The average front and back vertical wind 
shears in the direction of storm motion are 7, -- 1 
and -6 m s ^ for left-, straight- and right-moving 
cyclones, respectively. The same shear values at 
T - 36 are 3, 0 and - 10 m s'^ (see Fig. 14). Thus, 
for left-turning storms the net vertical wind shear in 
the direction of storm motion in octants I (front) 
and 5 (back) is positive. The shear is about zero 


for a straight-moving storm and negative for a right- 
turning storm. 

Such differences in tropospheric vertical wind 
shear between the turn classes only exist 24 and 36 h 
prior to turn time (Fig. 15). This shear difference 
does not show up as the storm approaches its turn 
time. Windshear values at T - 48 and before cannot 
be obtained because of the lack of cases. 

Thus, if one is able to obtain winds in the upper 
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FiCi. 14. Cross-sectional view of the sum of average V> wind components (m s M within 7- IP radii 
from the storm center in octants I and 5 at 200 and 900 mb at 24 h (T - 24) and 3b h (1* 3b) betore 

turn time for the three turn classes. The circled numbers are the values of the vertical wind shear 
[V/4200 mb)-V,.(900 mb)], 


200 mb 7 


^DIRECTION OF 
^STORM MOTION 
900 mb 4 


JuNii 1980 


CHAN, GRAY AND KIDD DR 


787 


and lower troposphere around a storm (for example, 
from satellite observations), it may be possible to 
better anticipate cyclone turning motion 24-36 h 
before it takes place. 

individual case analysis 

To see how valid the composite results are for 
individual turning storms, the Vp vertical shear de- 
scribed in the previous section for individual cases 
in the sample was aiso studied. Individual case 
analysis is severely handicapped by data avail- 
ability. Wind and height data must be interpolated 
from analyzed weather maps. Nevertheless, an at- 
tempt was made to evaluate this vertical shear 
parameter [V/j( 200 mb)-VV(900 mb)] for individual 
turning situations. The 200 mb winds were estimated 
from the height contours, using the geostrophic 
wind relationship. If there were no contour maps 
an estimate was made from the 200 mb streamline 
analysis. If neither was available, the case was not 
analyzed. The 900 mb winds were estimated from 
either the 3000 ft PIBAL maps or from the NHC 
ATOLL (top of the Ekman Layer) maps. If neither 
was available, an estimate was made from surface 
pressure maps, again using the geostrophic wind 
relationship. These methods restricted the analysis 
of the Vp vertical wind shear to only about two- 
thirds of the individual cases in the sample. 

Fig. 16 shows the shear values obtained at both 
T - 24 and T - 36. Although there is a large spread 
in the vertical wind shear for each classification 
(which is to be expected because of the inherent 
data uncertainties) a relationship between the verti- 



Fig. 15. Vertical wind shear for different time periods. The 
ordinate is the value of the average vertical wind shear (m s M 
within 7-11° radii from the storm center in octants 1 and 5. 
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Fig. 16. Vertical wind shear between 200 and 900 mb at T - 24 
and T-36 time periods for individual storms. The ordinate is 
the same as in Fig. 15. Each cross represents the wind shear 
value from one individual storm. Equal values of the wind shear 
are plotted next to each other. The dashed lines indicate the 
threshold values described in tb^ text. 


cal wind shear and the turning motion, similar to 
the one in the composite, does show up, 

To make sure that the results for these individual 
cases are not due to random chance, a / test was 
performed on the three samples of individual storms. 
The results for both T 24 and T-36 are shown 
in Table 3. The symbols are defined as follows: 

N number of cases in which the vertical wind 
shear parameter can be estimated 
X sample mean vertical wind shear 

jx average vertical wind shear of all three cate- 

gories of storms from composite data 
(assumed to be the population mean) 
s standard deviation of the sample 

/(959f) value of / at 95% significance level given in 
Panofsky and Brier (1958). 

It is seen that the mean vertical wind shears for 
both the left- and right-turning cases are significant 
at the 957c level. This lends confidence to the 
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Tabu. 3. The t test of the vertical wind shear parameter on the three samples of individual storms for both 
T 24 and T - 36. Symbols are defined in the text. 


Turn 

classification 

N 

X 


s 

t = |,v - mI<-V ~ 1)"= s ' 





T 24 




Left-turning 

U 

7.18 

0 

8.67 

2.62 

'*.23 

Straight-moving 

20 

-0.45 

0 

4.54 

0.43 

2.09 

Right-turning 

18 

-8.72 

0 

5.07 

7.09 

2.11 




T - 36 




Left-turning 

11 

6.23 

-2.33 

10.27 

2.64 

2.23 

Straight-moving 

14 

-L7I 

-2.33 

4.60 

0.49 

2.16 

Right-turning 

15 

“10.33 

“2.33 

6.39 

4.68 

2.15 


conclusion that tropospheric vertical wind shear can 
be used as an indicator of turning motion 24-36 h 
before such motion occurs. 

If we set threshold values for the 200-900 mb 
Vp vertical wind shear as 

>5, Left turn 
-5 to 5, Straight 
< “5, Right turn 

it is found that 64% of the left-turning storms have 
vertical wind shears > 5 for both T - 24 and 
T - 36 time periods. For the straight-moving 
storms, 77% occurred within the range at T - 36 
and 75% at T - 24, For the right-turning storms 87% 
have vertical wind shears < -5 at T - 36 and 72% 
at T 24. The present level of success in fore- 
casting turning motion is considered to be signifi- 
cantly less than this, as evident from the scatter of 
forecast points shown in Fig. 4. 

Since the techniques in deriving lower and upper 
tropospheric winds from (geostationary) satellite 
observations are improving, the information on 
such vertical wind shear should become more avail- 
able as time goes on and there may be hope for an 
improvement in future turning motion forecasts. 


6, Study of satellite temperature sounding data 


The vertical wind shear structure described in 
Section 5 should be related through the thermal 
wind relationship to the mean tropospheric tempera- 
ture differences across these storms. An estimate 
of the magnitude of the temperature gradient can be 
obtained from the thermal wind equation 


A/i R 





(1) 


where A/7 is the horizontal distance across which a 
temperature difference AT is measured, V p is the 
thermal wind, or in this case, V/.(200 mb) ~ V/,(900 
mb), /?, = 900 mb, = 200 mb, / is the Coriolis 
parameter and R the gas constant. If we assume / 
to be 5 X 10"*"' (corresponding to 20°N) and 


A/7 - 18® latitude, then for a vertical wind shear of 
10 m s ', AT ^ 2.2 K. 

To investigate the existence of such a horizontal 
temperature gradient, temperature sounding data 
from the Nimbus 6 Scanning Microwave Spectrom- 
eter (SCAMS) of 1975 were examined. Since this 
sounder can resolve tropospheric temperature gradi- 
ents of the above magnitude, we should expect the 
sounder data to be able to help specify turning 
motion. 

a. The Scanning Microwave Spectrometer 

The SCAMS on board the Nimbus 6 polar-orbit- 
ing satellite is a five-channel instrument sensing 
radiation nominally at 22.235, 31.65, 52.85, 53.85 
and 55.45 GHz (Staelin et a!., 1975a). The latter 
three frequencies are used for sounding the atmos- 
phere, while the first two are used to infer vertically 
integrated atmospheric water vapor and liquid water 
content. The Nimbus 6 spacecraft is in a sunsyn- 
chronous (1130 and 2230 Local Time) near-polar 
orbit at an altitude of approximately 1 100 km. The 
SCAMS instrument scans across the spacecraft 
track in 13 steps at increments of 7.2° in the scan 
angle each 16 s which results in a spatial resolution 
of 145 km near nadir, degrading to 200 km down- 
track by 360 km crosstrack at the maximum scan 
angle (Staelin et aL, 1975a). Approximately 400 
individual soundings can be obtained within 15° of 
the center of a tropical cyclone for each 12 h period. 

b. Tropical cyclones analyzed 

A total of 13 cases'^ (II individual storms) in both 
the Pacific (northeast and northwest) and Atlantic 
Oceans during 1975 were examined. As explained 
in Section 2a, only cyclones with maximum sus- 
tained wind ^ 18.0 m s ' were included. The tracks 
for these storms are shown in Fig. 17. The distribu- 
tion of storms in each turn class is given in Table 4. 
Satellite data were categorized by time (T, T ~ 12, 


This is all the SCAMS information available to the authors 
at this time. 
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'Faui I 4. Dislrilnuion of turn discs in the Atlantic and Pacilic 
Oceans used in studying satellite temperature sounding data 
around tropical cyclones. 


Ocean 

l.eft 

Straight 

Right 

Atlantic (West Indies) 

0 


4 

Northwest Pacilic 

1 

I 

3 

Northeast Pacilic 

1 

1 

0 

Total 

2 

4 

7 


T -- 24 and T - 36) as described in Section 3. The 
mean 1000-250 mb temperatures around each 
storm for each time period were derived from the 
brightness temperatures (the raw data obtained from 
the sounder) as described in Waters ct al, (1975). 
'Fhe same method used in the rawinsonde composite 
was employed (see Section 3) to obtain average layer 
mean temperature at each of the grid points in J^'ig. 7 
for each storm. To minimize contamination by 
precipitation, data points with integrated liquid 
water contents 0,5 kg m - were not used in the 
averaging process. According to Staelin ct cil. ( 1975b) 
this should eliminate most of the contaminated data. 

The accuracy of atmospheric temperatures de- 
rived from satellite sounders is a pertinent point. 
Data from a microwave sounder were used because 
microwave soundings are unaffected by non-raining 
clouds (Staelin ct aL, 1975b). The extensive cloudi- 
ness near tropical cyclones precludes the exclusive 
use of infrared sounding data. When the 1000- 
250 mb layer mean temperature at a single spot is 
compared with the nearest grid point data analyzed 
by the National Meteorological Center (NMC) the 
root-mean-square difference is --2 K (LandsaV 
Nimbus Project, 1976). It must be rcmcmhcrccL 
however, that only temperature difference and not 
absolute temperature is a requirement for the wind 
shear determination. In addition, each grid point 


value is an average of a considerable number of 
observation points, and precipitation-contaminated 
data points have been eliminated. The satellite- 
derived temperature differences obtained in the 
way described above should thei’efore be accurate to 
within a few tenths kelvin (Staelin et aL, 1975b). 

c, 1 000 -‘250 mb mean temperature Jields 

For convenience, the l()()0-250 mb layer mean 
temperatures were analyzed in the following way. 
All layer-mean temperatures at grid points within 
5-15' radius’’ from the storm center were averaged 
to obtain an average layer-mean temperature 1\ 
The deviation of each grid point temperature from 
t was calculated (for all points in the grid), These 
deviations were used for the analysis instead of the 
absolute values. An example of the temperature 
deviation field for three individual cases 24 h 
before turn time (one in each turn class) is shown in 
F\(^, 18. This figure shows the typical expected differ- 
ence in temperature structure existing between a 
left-turning, straight-moving and right-turning storm 
as implied by the difference in vertical wind shear 
structure discussed in vSection 5. The left-turning 
storm shows a temperature gradient across the 
storm corresponding to a thermal wind, or vertical 
wind shear, in the direction of storm motion while 
the opposite occurs for the right-turning storm. 
There is very little temperature gradient for the 
straight-moving case. An examination of the other 
cases indicates that the temperature gradients of 
most significance exist in front of the storm. In 
addition, for most of the left- and right-turn cases 
studied, there exists a region where the troposphere 


The reasons for using 5-15^" are !) more data are available 
outside 5‘’ (a lot of data were eliminated within 5 ' because of high 
liquid watercontent), and 2) 15’is the outermost radius used in the 
rawinsonde composite study. 



Fio. 18. 1000 to 250 mb layer mean temperature deviation from 5 to 15^ radii average mean temperature derived from the Nimbus 6 
SCAMS temperature sounding data 24 h before turn lime. The arrow from the sUirm center indicates the direction of storm motion. 
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T,\in i: 5. The 1000-250 mb layer mean temperature Oevialions from 5-- 15' radii average in octants 2 (len front) and 8 (right front) 
within 9- UT‘ radius from the storm center for various storms tn the North Racilic and Atlantic Oceans 24 h before turn time. Vertical 
wind shear in the direction in which the storm is going is calculated from the thermal wind relationship using a horizontal distance 
of !8\ The Coriolis force./ is taken to be 5 10 s 


Storm 

Region 

Octant 2 

Octant K 

Octant 2 
minus 
Octant 8 

Corresponding 
vertical 
wind shear 
(m s *) 

Turn 

Phvllis 

west Pacilic 

-0.8 

0.2 

1.0 

4.5 

left 

lisa 

east Pacilic 

0.9 

1.5 

2.4 

10.9 

left 

Alice 

west Pacilic 

0.7 

0.0 

0.7 

3.2 

straight 

Lily 

east Pacilic 

- 0.6 

0.1 

0.5 

2.3 

straight 

Kloise 

west Atlantic 

0.2 

0.2 

0.0 

0.0 

straight 

Caroline 

west Atlantic 

-LI 

2.9 

-4.0 

18.2 

straight 

Hloise 

west Atlantic 

1.5 

0.2 

1.7 

-7.7 

right 

Hloise 

west Atlantic 

1.6 

-2.0 

3.6 

-16.4 

right 

Gladys 

west Atlantic 

0.1 

-1.4 

1.5 

6.8 

right 

Faye 

west Atlantic 

0.3 

0.4 

0.7 

-3.2 

right 

June 

west Pacilic 

0.7 

-1.6 

2.3 

- 10.5 

right 

Tess 

west Pacilic 

-0.7 

- 1.9 

1.2 

5.5 

right 

Cora 

west Pacilic 

1.0 

1.0 

0.0 

0.0 

right 


is relatively cold in the direction toward which the 
storm is moving, as indicated in the examples given 
in Fig. 18. 

To describe these temperature gradients in more 
quantitative terms, temperature deviations at 
9~I3° radius*’ in octants 2 (left front) and 8 (right 
front) for all cases were studied. Results are shown 
in Table 5. Both left-turning cases show that the 
mean tropospheric temperature is lower to the left 
front of the storm. For the straight-moving storms, 
only Caroline shows a large negative tropospheric 
shear. The reason for this is unknown. However, all 
other straight-moving cases show only a very weak 
temperature gradient across the front of the storm. 
In addition, all three cases show a region of cold air 
in octant 1 (not shown). Of the seven right-turning 
storms, only Cora shows no temperature gradient. 
The others indicate a relatively cold troposphere to 
the right front of the storm. 

These results were obtained, as described earlier, 
after elimination of data points with integrated 
liquid water content 0.5 kg m To see if this 
process is of any great significance, temperature 
deviation fields without eliminating any data were 
also analyzed. A comparison between the two sets 
of results shows very little difference in the tempera- 
ture structure and temperature gradients at large 
radii from the storm center. It would appear that it 
is not necessary to eliminate those data points 
with appreciable liquid water content provided 
temperatures at large distances from the storm 
center are used. 


In order to obtain the horizontal temperature gradient normal 
to VV in octant I, at radii 1 - 1 1 (the area where the vertical wind 
shear parameter is of signilicance), one has to consider tempera- 
tures in octants 2 and 8 at radii 9-13'’ (see Fig. 7). 


These results suggest that very pertinent informa- 
tion on forecasting tropical cyclone turning motion 
24 h in advance may be obtained from satellite 
sounder information in region*^ devoid of rawinsonde, 
aircraft and clouds for satei*'i'. wind vector deter- 
mination. Newer and better satellite sounders are 
likely to be available in the future for even better 
determinations. 

(). Summary and recommendations 

The results of this study indicate that there are 
parameters in the large-scale environment around 
tropical cyclones at upper and lower tropospheric 
levels (in addition to the previously considered 
crucial middle level) which can provide infoimiation 
about their turning motion 24-36 h before it occurs, 
An estimate of the V'y component of the winds at 
200 and 900 mb at 5- IF radius around the storm 
from satellite observations may also be useful for 
making short-range (0-12 h) forecasts of turning 
motion. Of greater significance is the possibility of 
using the F/- vertical wind shear to specify turning 
motion 24-36 h in advance. T'he fact that favorable 
environmental conditions are evident 24-36 h prior 
to turn time indicates that there is likely a time lag 
between changes in the environmental fiow param- 
eters surrounding a storm and the response of the 
storm to such changes. 

At certain times it may not be possible to derive 
satellite winds at 900 and 200 mb because of an 
extensive cirrus canopy or because upper-level 
cloud motions are lacking. In these cases the satel- 
lite sounder, with its tropospheric temperature 
difference measurement capability, may be useful in 
augmenting missing wind data. vSince contamination 
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of the sounding data by precipitation ivS negligible 
at large radii, the Microwave Sounding Unit (which 
does not have any water vapor channel) on board 
the TIROS-N satellite now in orbit might be used to 
monitor tropospheric temperature gradients across a 
storm to determine the upper and lower tropospheric 
vertical wind shear. It is suggested that temperature 
data obtained from such a sounder be used on an 
operational basis to test the applicability of these 
results. 

With the large increase in computer technology 
and the development of more sophisticated cyclone 
forecast schemes, one would expect a continuous 
increase in tropical cyclone forecast skill. While 
this has been true on the average in the 1960\s and 
early 1970\s, tropica! cyclone track 24 h forecast 
errors have been showing a slight increase in the 
last three to four years (see, e.g,. Joint Typhoon 
Warning Center, 1977). These forecast error in- 
creases are likely a result of an over-reliance by the 
various forecast techniques on middle-level flow 
patterns, and a subsequent gradual reduction of 
oceanic rawinsonde and aircraft data to measure 
such flow patterns. Atlantic weather ships Echo 
and Delta were decommissioned in 1973. Very few 
oceanic propeller aircraft flights are now made. 
These changes have significantly reduced the 
already scanty mid-tropospheric Atlantic wind- 
height data. Accurate analysis of middle-level flow 
fields has always been difficult. At some locations 
and time periods it is now nearly impossible to 
obtain a meaningful analysis. Neumann (private 
communication) has discussed how NMC tropical 
analysis schemes can produce a quite unrepresenta- 
tive middle-level steering current determination. 
He has documented cases in which an erroneous 
wind observation drastically changes middle-level 
computer-analyzed flow fields. This seriously 
affects the performance of tropical cyclone forecast 
models based on middle-level steering flow concepts 
(e.g., see Lawrence, 1979). New forecasting schemes 
(regardless of how sophisticated) which utilize 
middle-level flow patterns should encounter similar 
observational problems. 

Because of the increasing availability of satellite- 
derived information at upper and lower tropospheric 
levels, it is suggested that research on storm mo- 
tion be concentrated on studying the relationships 
between upper and lower tropospheric parameters 
and storm movement. The next steps in this research 


will be to try to relate upper and lower level 
parameters with cyclone speed changes, stalling and 
looping motions, and other special cyclone motion 
characteristics. 
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ABSTMCT OF THKSIS 


SUMMER PRECIPITATION FREQUENCY OF THE NORTH CENTRAL U.S. FROM 
SATELLITE MICROWAVE OBSERVATIONS 

The Niinbus-6 Electrically Scanning", Microwave Radiometer (ESMR--6) 
is used to observe raining areas over the north central U.S. from 11 
June through 31 August 1976. A technique is developed whereby the 
frequency distribution of microwave brightness temperatures over the 
area is analyzed to determine a threshold temperature below which 
the possibility of precipitation is assumed. Raining areas are clas- 
sified by considering spatial distribution of possible raining areas, 
simultaneous infrared Images, surface temperatures, surface relative 
humidities as \<^ell as microwave brightness temperatures. Preci])itation 
frequency of the area derived by the satellite technique is compared 
to surface raingage hourly observations. The relationship between the 
two systems of measurement is seen to be quite high in the southern 
and northwestern quadrants of the area however the northeastern quadrant 
exhibits poor relationship due to the existence of many small lakes in 
Minnesota and Wisconsin. Precipitation frequency determination from 
satellite microwave observations is feasible for land areas warm enough 
to be distinguishable from raining areas. Thus, determinations over 
primary mid latitude crop growing seasons are possible. 

Philip A. Durkee 
Atmospheric Science Department 
Colorado State University 
Fort Collins, Colorado 80523 
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1.0 INTRODUCTION 


Rainfall measurement is important to many disciplines i A^^riculturc 
needs to know precipitation patterns and amounts to determine irrigation 
needs and flooding potential. Hydrology and Water Resource Management 
needs rainfall measurements to determine soil moisture, runoff potential 
and reservoir allotment. Knowledge of rainfall is important to the 
atmospheric sciences in understanding latent heat release in the at- 
mosphere as well as determining climate tendencies. 

The measurement of rainfall over land areas has primarily been 
done by raingages although radar observation of precipitation has 
been an important technique. Spatial coverage by these measurement 
techniques is limited, however. A network of raingages observes single 
point values of the precipitation over an area and interpolation be- 
tween gages is highly approximate. Radar observations, though over 
a continuous area, have only a finite range to which the instrument can 
reach. Also, both observing systems have inherent errors since arrays 
of different instruments each with different operators are assumed the 
same in analyses. 

The development of satellite techniques has aided greatly in 
eliminating some of the problems inherent in surface observation of 
precipitation. Griffith, et al. (1978) observed the growth of convec- 
tive clouds using geosynchronous infrared and visible imagery to esti- 
mate rainfall. The use of microwave frequencies allows direct obser- 
vation of precipitation due to the large interaction of microwaves with 
hydrometers (falling raindrops) and small interaction with clouds and 
gaseous constituents (away from ^60 GHz). The Nimbus 5 Electrically 
Scanning Microwave Radiometer (ESMR-5) , operating at 19.35 GHz, has 
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been used to observe oceanic precipitation ^'Kidder, 1976; Adler and 
Rodgers, 1977; and Wilheit, et al. , 1977), an area with little or no 
surface observations. Rainfall observation over land areas was shown 
feasible by Savage and Weinman (1975) using the Nimbus 6 ES^^R operating 
at 37 GHzo Investigation of ESMR-6 rainfall estimation was continued 
by Savage (1976), Weinman and Guetter (1977), Weinman and Davies (1978), 
and Rodgers, et al. (1978). Precipitation frequency of the north cen- 
tral U.Sc in summer using ESMR-6 observations is the subject of this 
paper. 

The theory of microwave radiative transfer is presented in Chapter 
2 including a solution of the radiative transfer equation, absorption 
by water vapor, liquid water and oxygen, and a consideration of scat- 
tering by hydrometers. Chapter 3 describes the ESMR-6 instrument and 
the data sets used in this study. Also included in Chapter 3 is a 
description of the thresholding technique used to classify raining 
areas within the area of interest. Surface emission and the relation- 
ship between the satellite observed brightness temperature and surface 
temperature are also examined. Chapter 4 presents the results of the 
precipitation analysis as well as a comparison with hourly raingage 


observations. 


2.0 MICROWAVE RADIATIVE TRANSFER 

2.1 The Radiative Transfer Equation 

For a plane-parallel atmosphere the equatioii of radiative transfer 
may be written (Chandrasekhar, 1960) 
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dlv 

dz 


■ -0 


(1) 


where y = cosine of the zenith angle 

" radiance at frequency v at position x at time t 

z = vertical coordinate 

p = mass density of the atmosphere at x 

Kv = volume absorption coefficient per unit mass 

= S (x,t) = source function of a unit mass at position 
X at time t 

00 

It \-L pKvdz’ is the normal optical depth then (1) can be written 
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Multiplying (2) by exp( — jj-) gives 
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The left hand side of (3) is just 
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so integrating (3) with respect to optical depth from the surface to 
Ty we get 


I 

V 





(5) 


is the radiance at frequency V out the top of the atmosphere which 
may be viewed by a satellite. 

The source function, S^, is a function of the radiative properties 
of the atmosphere. For a scattering, non-absorbing atmosphere 

\ Iv sineMeMcj)’ (6) 
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where 0 and ^ are the zenith and azimuth angles respectively. p(0,0’,Ct^*) 
is known as the phase function and represents the angular distribution 
of scattered energy. For a non-scattering atmosphere satisfying local 
thermodynamic equilibrium 
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where h = Planck’s constant 

k = Boltzman’s constant 

c = velocity of light 

T = local thermodynamic temperature 

h / 

and is just Planck’s black-body relation. If — <<1 then e may 

be expanded and may be approximated as 

2v\t 
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( 8 ) 


5 


which is the R.ayleigh-Jeans approximation. At 37 GHz (which is the* 
frequency of interest for this study) for example hv/VT = 0.006 if 
T = 300 K v?hich introduces an error of less than 0.3% (Kreiss, 1968). 
Defining the brightness temperature, such that 


I . = T, 


(9) 


the solution of the radiative transfer equation becomes 
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T, = e ^ + 
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( 10 ) 


The contribution to the total brightness temperature by the 
surface, is made up of two parts; the emitted radiation and re- 

flected radiation incident downward on the surface. The satellite ob- 
served brightness temperature from equation (10) is then just the portion 
of the effective surface brightness temperature which is transmitted 
through the atmosphere, plus the contribution from atmospheric emission 
which is transmitted. 


2.2 Transfer of Microwaves Through the Cloudy Atmosphere 

The three major atmospheric constituents which absorb and hence 
emit thermal microwaves are molecular oxygen, water vapor, and liquid 
water. Figure 1 shows the absorption coefficient near the surface of the 
earth of oxygen and water vapor as a function of frequency. Water vapor 
absorption in this region is due to a rotational transition at about 22.2 
GHz. Oxygen absorption results from a series of magnetic dipole moment 
vs. rotational angular moment reconfigurations from 53 GHz to 66 GHz 
(Meeks and Lilley, 1963). Figure 2 shows the absorption coefficient of 
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MICROWAVE ABSORPTION COEFFICIENTS AT STANDARD 
SURFACE CONDITIONS 


Figure 1 



Absorbtion coefficient for atmospheric water vapor 
and oxygen. (After Westwater, 1965) 
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various liquid water contents of the atmosphere as a function of 
frequency from Staelin (1966). The relationship is valid from 10 
to 33 GHz and is a good approximation for small droplets compared to 
X/2 'jt, where A is the wavelength. As droplets increase in size, their 
scattering effects become more important. These eifects will be 
discussed in section 2.4. 

Of the microwave absorbing constituents of the atmosphere, water- 
vapor and liquid water are the most variable. These also have a sig- 
nificant, though small, temperature dependence on absorption. Figure 
3 illustrates this temperature dependence for various mass densities of 
water vapor and liquid water. These curves are calculated by expres- 
sions from Van Vleck (1947) and Staelin (1966) for water vapor and 
liquid water respectively. Both show a decreasing absorption coef- 
ficient with increasing temperature. 

The temperature dependence of x^rater vapor absorption is due to 
the temperature dependence of the Boltzman distribution. Simply it 
is a result of the fact that with higher temperature fex^rer molecules 
will have an energy corresponding to the lower energy state of the 
transition resulting in less interaction xijrith the radiation and less 
absorption. As apparent from Figure 3 the temperature dependence is 
small for water vapor absorption. For a decrease of 40 K in tempera- 
ture, the absorption coefficient increases by 30%. Temperature depen- 
dence on liquid water absorption is due to the temperature dependence 
on the dielectric constant of water. As temperature increases the 
dielectric constant of x^ater decreases « The result from Figure 3 is 
that with a 40 K decrease in temperature the absorption coefficient 
increases by a factor of 3. 
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It should bo noted that xoe particlor* are hi^Uily transparent at 
microwave frequcmcieH, Water covered ice however, has ln»en sliown to 
have significant absorbtive properties (Westwater, 1972)* 

The total optical depth, is obtained by integrating the 

-- Iv 

absorption coefficient with height. Tlie factor e h in (3) is then 
kno\fn as the fractional transmittance of the atmosphere. Calculations 
by Westwater (1972) show that at 37 the clear atmosphere has a 
transmittance of about 93% for the U.S. Standard atmosphere and about 
88% for the 15^N Tropical atmosphere. Calculations were also made 
using typical cloud liquid water and water vapor distributions for 
non-precipitating clouds. Trausmittances remained greater than 70%. 


2.3 Surface Emission 

At microwave frequencies below 50 GHz the atmosphere is liighly 
transparent. Therefore the effective brightness temperature of the 
surface, ^^^9 often the most prominent contribution to the brightnes 
temperature observed at the satellite: 


bo 


- 6 T^, 
V s 


+ T. 


hr 


( 11 ) 


where € and T are the emittance and temperature of the surface and 
V s 

T, is the contribution due to reflection of incident radiation. The 
br 

surface emittance then is the important characteristic which controls 
a large part of the total signal. 

In the microwave region surface emittances are a function of : 
terrain composition, electrical conductivity of the surface, roughness 
of the surface, and the frequency of observation. At 37 GHz the 
emittance of a water surface is about 0,45 while for most land sur- 


faces the emittance is about 0.90. 
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Since the transmittance of tlie earth’s surface is zi^roy tlu* 
omittance is related to the refiectivt» properties of a surface hy 

e - I - r , ( K! ) 

where v is the frequency and r^ is the reflectance of the surface. 

From (12) surfaces of high omittance have low reflectance and low 
emittance infers high reflectance. It is the relatively high re- 
flectance of water surfaces which distinguishes them from raining 
areas over land surfaces. 

Radiation incident on a reflecting surface at some angle will 
be partially polarized to the extent that the angle of incidence is 
close to Brewster’s angle for the medium. For water surfaces Brewster’s 
angle is very close to 50^ which is the constant viewing angle of 
the ESMR-6 , Therefore brightness temperatures from water surfaces will 
be lower than from land surfaces (due to lower emittance) and will also 
be relatively polarized (due to higher reflectance) with vertical 
polarization greater than the horizontal polarization. 

2.4 Scattering of Microwaves by Hydrometers Over Land Surfaces 

As cloud droplets increase, in size and approach precipitation 
sized particles the effect of scattering of microwave radiation in- 
creases. Deirmendjian (1963) has showni that at frequencies below 60 
GHz, for clouds with a modal drop size of greater than 50 \m scat- 
tering effects become important compared to absorption. In a scattering 
atmosphere equation (G) must be used in addition to the Planck relation 
as the source function in the radiative transfer equation along witli 
the full Mie equations. Weinman and Guetter (1977) have calculated 37 
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GHz microwave radiances from rain clouds over land and water suriaces 
in both vertical and horizontal polarizations. Fij^ure 4 is a plot of 
their results for various rain rates over land, roup,h" and calm water. 
Brightness temperatures of rain clouds is seen to be a decreasing 
function of rain rate over land surfaces x^hile over water surfaces 
brightness temperatures increase with rainfall rate to about 4 mm/hr 
and then decrease as over land. Above 8 mm/hr the surface is no longer 
an important factor in determining the brightness temperature as the 
emission of the raining cloud is the sole contributor to the upward 
radiance. 




The distinction of a rough water surface is only in that the 
surface does not polarize. 


BRIGHTNESS TEMPERATURE (°K) 


CALM WATER 

ROUGH WATER 

DRY LAND 



Figure 4. Brightness temperatures of 37 GHz radiances as a 
function of rainfall rate over land, rough and 
calm water. (After Weinman and Guetter, 1977) 


3.0 DATA AND PROCEDURES 


3.1 Description of the Instrument 

The Nimbus 6 satellite was launched on 12 June 3975 into a 
circular sun-synchronous polar orbit at an altitude of approximately 
1100 km. On board was the Electrically Scanning’; Microwave Radiometer 
(ESMR-6) which detected the vertical and horizontal polarizations of 
thermal microwaves in a 250 MHz band centered at 37 GHz. ESMR-6 
operated satisfactorily through 15 September 1976 when the horizontal 
channel failed* The instrument was situated aboard the satellite as 
shoxm in Figure 5 where it electrically scans along a cone forward of 
the satellite such that the incidence angle at the surface remains 
constant at 50°. The scan is done in 71 discrete steps from 35° to 
the right to 35^ to the left of the forward direction. For a zero 
degree pitch angle the resolution of each beam position is 20 Ion 
cross-track and 42 Ion down-track. 

The orbit of Nimbus 6 had a retrogade inclination of 81° with 
equator crossings at local noon (northward) and midnight (southward) 
and a period of 107.25 minutes. Figure 6 shows a typical 24 hour 
coverage of the ESMR-6. The width of the 71 beam position scan is 
only 1272 km which results in large gaps in coverage over the equa- 
torial regions. The gray shaded areas and areas between 60° and 85° 
latitude are covered at least twice each 24 hour period. There is 

O it 

no coverage polex^ard of 85 . 


See Wilheit (1975) for a complete description of the ESMR-6. 
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Figure 5. The Nimbus 6 satellite. (From The Nimbus 6 User’s 
Guide) 
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3.2 ThG ESMR-6 Data Set 

The ESMR-6 data used in tliivS study was recorded from 11 June 1976 
through 30 August 1976 over the north central U.S. as shomi in Figure 7 
Prior to this period the ESMR was operated only on a half-time basis 
while after 15 September 1976 the horizontal polarization was unavail- 
able. This period then represents the most intensive temporal coverage 
of the ESMR-6. Figure 8 shows the distribution of available orbits 
over the entire area, so that any one point in the area is viewed by 
the satellite only about 35 times out of 57 noon orbits and 40 times 
out of 72 midnight orbits. 

3.2.1 Corrections Made to this Data 

The microwave brightness temperatures are recorded on magnetic 
tape and available from the National Space Science Data Center (NSSDC) 
at the Goddard Space Flight Center (GSFC) . Processing of the data at 
NSSDC includes calibrations which correct for thermal effects of the 
instrument, sidelobe effects, and polarization mixing. To investigate 
the completeness of this calibration procedure the data to be used in 
this study was averaged by beam position for vertical and horizontal 
polarization of noon and midnight orbits. The results of this averaging 
are shown in Figure 9. If all necessary corrections had been taken 
into account, the brightness temperature would be the same for each 
beam position in a long term average. As is seen in Figure 9 large 
variations of average brightness temperature with beam position exist, 
particularly in the outer 5 beam positions. These large variations 
in the outer positions are due to interference from the sidelobes of 
the phased array antenna. Because of the large and uncertain magnitude 
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Figure 7» Areal coverage of this study. 
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Figure 8. Distribution of available ESMli-6 observations used 
in this study. 
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Figure 9 




Average brightness temperature for each beam position 
along the ESMR-6 scan over all observations used in 
this study, (a) noon observations, (b) midnight 
observations. 
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of thu Kidtilobe beam posiLiouH I-'; and 07-71 aru ukgIuss for 

quantitative analysis. The variations over the rcnaainin^^ beam positions 
are probably due to imperfect polarii'.ation mlxiny, corrections and slipjit 
lobe interference. To account for these variations an empirical correc- 
tion was applied to all beam positions according to the difference of 
the beam position's average brightness temperature from the average 
brightness temperature of the center beam position (position 36) . Table 
1 lists the values subtracted from the brightness temperatures by beam 
position, polarization and time. 


3.3 The Surface Raingage Data Set 

In this study hourly rainfall observations are used for comparison 
with the satellite derived rainfall. These observations were obtained 
from Hourly Precipitation Data . The stations for which surface raingage 
observations are available are sho\>m in Figure 1C. 

Rainfall occurrence has been tabulated for the hour of each pass 
of the satellite for those stations which lay within the area of ESMR-6 
coverage. The results of the surface raingage observations are given 
in Chapter 4. 


3 . 4 Thresholding Technique 

' For a non-scattering atmosphere, the microwave brightness tempera- 
ture observed from a satellite may be writteii (Westwater, 1972); 


■/ 


>c - fo 

m •'O 

Tae 


’ds’ 


ds 


+ € T 
s s 




(13) 


where t = 



CO 

a ds is the total optical depth, i and s are displacements 


Beam Posi 
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upward and downward ruspoctivaly, 6 and T art* the omittanoo anci 
temperature of the surface, T is the temperature as a function of 
height and y = cost) where 0 is the zenith angle. then is simply 

the sum of contributions from the emission of the atmosphere between 
the satellite and the surface, the emission of the surface and the re- 
flection of the amission of the atmosphere incident downward on the 
surface. If we assume that total brightness temperature is a function 
of surface emittance and temperature only, that is variations in the 
vertical temperature profile and distribution of absorbing and emitting 
constituents are negligible, then equation (11) may be written 

- A -f e T B + (1 - € )C (14) 

where A, B and C are constants. Westwater (1972) has made calculations 
of the contributions to T, for both the U.S. Standard Atmosphere and 

D 

the 15^N Tropical Atmosphere. Table 2 lists the interpolated values 
of the terms in equation (14) for = .9 and - 288 for the U.S. 
Standard Atmv ^^here and T^ = 299 for the 15^N Tropical Atmosphere. It 
is apparent from Table 2 that at 37 GHz the dominant term in deter- 
mining the total brightness temperature is the surface emission which 
accounts for 88% of the U.S. Standard Atmosphere total. Upward emission 
of the UoS. Standard atmosphere accounts for 12% while only 1% of the 
total brightness temperature is contributed by the reflected downward 
emission of the U.S- Standard Atmosphere. The 15^N Tropical Atmosphere 
is more moist and warmer than the U.So Standard Atmosphere whicli results 
in a larger portion of the total brightness temperature being contrib- 
uted by the atmospheric term. The result is that 18% of the total is 
added by upward atmospheric emission, 80% is added by the surface 


Numerical values of the tormr. in Kciuation 
from Westwatcr (1072) 

(14) 

U.S. Standard 

1')"n Tropical 

30.0 K 

49.5 K 

0.833 

0.823 

30.0 K 

49.5 K 

229.0 K 

221.5 K 

3.0 K 

5.0 K 

262.0 K 

276.0 K 
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emission and 2% of the total is added by the reflected downward emission 
of the 15^N Tropical Atmosphere, The results of this partitioning, of 
the total brightness temperature show that the variations in the surface 
terms, emittance and temperature, will account for the ma;jority of the 
variance in satellite observed brightness temperatures over land if 
large areas of scattering hydrometers are not present. 

From Table 2 (term C) it is seen that the brightness temperature 
of the downward radiation ixicident on the surface is about 30 K and 
49.5 K for a U.S. Standard and IS^^N Tropical Atmosphex.e respectively. 

If water surfaces have an emittance of 0.5 then the additJ of this 
term as seen by the satellite is about 15 K for a U.S. Standard atmo- 
sphere and 25 K for a 15^N Tropical Atmosphere. Nearly all of this ad- 
ded brightness temperature will be vertically polarized so that the 
brightness temperature in the vertically polarized channel will be 
higher than in the horizontally polarized channel. Therefore, to elim- 
inate cold brightness temperature water surfaces from consideration, 
areas with greater than a 15 K difference between the vertical and 
horizontal channels are eliminated throughout this study. 

Figure 11 shows the frequency distribution of observed microwave 
brightness temperatures from the vertically polarized channel over the 
area in Figure 8 at 0630 UT on 14 June 1976. The shape of the distri- 
bution is that of a Gaussian with a tail extending from the cold end 
of the distribution. This distribution which is typical for this 
at this time of year, is comprised of two parts. First, the non-raining 
fields of view over the area make up the Gaussian part of the distri- 
bution with a characteristic standard deviation about some mean due to 
variations in atmospheric and surface radiative conditions. Second, 


NUMBER OF OBSERVATIONS 


14 JUNE 1976 0630 GMT 



TEMPERATURE (°K) 


Figure 11, Frequency distribution of observed microwave bright 
ness temperatures from the vertically polarized 
channel at 0630 GMT on 14 June 1976. 
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raining areas which decrease the observed brightness temperature of a 
field of view over land cause a displacement of temperatures towards 
the colder end of the distribution. The general Gaussian shape of the 
distribution is therefore distorted soniewliat on the cold end by fields 
of view partially to completely filled with raining elements. 

Rodgers et al. (1978) has shown statistically that indeed popu- 
lations of raining and non-raining fields of view are separable. If 
the standard deviation of the non-raining distribution is known, the 
distribution may be subtracted from the whole leaving the raining fields 
of view. 

The following numerical determination of the standard deviation 
and mean value of a contaminated Gaussian has been presented by Smith 
(1979). The functional form of a Gaussian distribution can be written 


£(i) . f )-Ct3£ 

max ^ 1 2 o *2 


xp 


(15) 


where f is the number of observations of temperature T, fmax the 
maximum number of observations for any one temperature, T is the 
average temperature or the temperature of maximum observation, and a 
Is the standard deviation. Given three points along the distribution 
say (f^, ip, (fy T^) and (f^^, such that then 


^ -1 
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which are two equations with the unknowns a and T. If a = 
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ai\d 
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2 (1-a) - Tj + uT,^ 


0 - 



(l-B) T,“ - t/ + UT. ^ 
1 j K 

(1.-3) Inf - Infj + 31nfj^ 


(lf<) 


(19) 


Equations (18) and (19) are applied to all possible combinations of 
three points which are safely in the non-raining portion of the dis- 
tribution, naraely the warm end of the Gaussian. The most common solu- 
tions then provide good estimates of a and T. The values of a and T 
derived for the 14 June 1976 distribution are indicated in Figure 11. 

The standard deviation is a measure of the width of the distribu- 
tion. It is therefore feasible to use some multiple of o halow the 
mean as the lower temperature boundary of the distribution of non- 
raining fields of view. Figure 12 is a plot of percent of the total 
Gaussian distribution below the mean minus some multiple of the stan- 
dard deviation vs. the multiple of the standard deviation. For a 
threshold of two standard deviations below the mean of the distribu- 
tion less than 2% of the total non-raining fields of view remain below 
the threshold. Since these Z'Z (on the order of 25 fields of view) are 
expected to have little spatial organization they will be eliminated 
from the climatology by only accepting connected fields of view below 
the threshold as raining areas. Also, misclassif ication of non-raining 
fields of view can be avoided by examining observed surface temperatures. 
If large areas of cold surface temperatures exist in the area, say 
behind a cold front,* fields of view below the threshold in these areas 
may be eliminated. 


Certainly along the cold front where precipitation is highly possible 
all connected fields of view would be accepted as raining areas in 
the climatology. 



Ivigure 12. Fraction of a Gaussian distribution less than 
imiltiple of the standard deviation of the dis 
bution below the mean of the distribution. 
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3,5 Surface Temperature vs. Brightness Temperature 

Frequency distributions of sur:ace temperature from the 43 first- 
order observing stations in the area ( Local Cli matological Data ) liave 
been compared to the satellite brightness temperature distributions to 
examine the relationship between the two quantities. Figure 13 is a 
plot of modal brightness temperature vs. modal surface temperature for 
the area of the study. Daytime cases are open circles and nighttime 
cases are closed circles and only orbits which covered more than 30 of 
the 43 surface temperature stations are shown. Assuming an omittance of 
0.9 yields the calculated curves for the 15°N Tropical and U.S. Standard 
atmospheric temperature and moisture profiles from Westwater (1972), 

The relationship for daytime cases appears to be well approximated 
by the 15°N Tropical Atmosphere. The nighttime carvm lie along the 
curve describing nearly surface emission. This infers that the night- 
time radiative conditions are such that the addition of atmospheric 
emission to the brightness temperature is in general cancelled by a 
decrease of brightness temperature due to the opacity of the atmosphere. 
The above conclusions are based on the assumption that oh • ^^-^ed surface 
temperature (which is actually air temperature) closely : i ates 

the average surface emitting temperature within a field Oj. " i •• )r 
that at least on the average over the area this is true. One would 
expect that daytime surface emitting temperatures would be somewhat 
higher than air temperature and that nighttime surface emitting tempera- 
tures would be somewhat cooler than air temperatures. If this were 
the case then the daytime and nighttime cases would be more co-linear 
along a line lying between the 15^N Tropical curve and the surface 


emisrion curve. 


MODAL SATELLITE BRIGHTNESS TEMP (°K 



MODAL SURFACE TEMP (°K) 


Figure 13. Modal brightness temperature vs. modal surface 

temperature. Open circles are daytime cases, closed 
circles are nighttime cases. (Only orbits which 
covered more than 30 of the 43 surface temperature 
stations are sho^smO 
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3.6 Classification of Raining Areas 

In preparing the precipitation frequency from the ESMR-6 ohsorvatiorib 
the following information sources were used: 

1) Microwave brightness temperatures: Areas where the brightness 

temperatures were less than 2 standard deviations below the 
mean of the distribution of non-raining brightness temperatures 
were taken as possible raining areas. 

2) Spatial distribution of possible raining areas: Only areas 

with 2 or more adjoining fields of view were taken as possible 
raining areas. 

3) Simultaneous THIR images: Images by the Temperature-Humidity 

Infrared Radiometer (THIR) (McColloch, 1975) also on board 
Nimbus 6 were used to assure the presence of clouds over the 
raining areas. THIR images were available for about 60% of 
the daytime ESMR-6 observations over the area of this study, 
while no nighttime THIR images were available. 

4) Surface temperatures: Surface temperatures observed at noon 

and midnight at the 43 first-order stations were examined so 
that large areas with cold surface temperatures could be dis- 
counted as raining areas. Identification of these areas were 
aided in that the fields of view with brithtness temperature 
less than the threshold were scattered throughout the area 
and showed little or no organization. 

5) Surface Relative Humidity: Midnight observations of surface 

relative humidity were examined so that areas with high proba- 
bility of dew formation could be discountea as raining areas. 
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Relative hunddities >90Z were considered to indicate the 
possibility of dew formation* Again, identification of these 
areas was aided by the unorganized spatial distribu+'ion of 
fields of view below the threshold in these areas. 

6) Polarization of the microwave brightness temperatures: As 

throughout this study only fields of view with low polarization 
(vertical channel not greater than 15 K higher than the hori- 
zontal channel) were considered. This eliminates water sur- 
faces from consideration. 

Figure 14 is a flow diagram illustrating the decision process used to 
classify fields of view as raining or non-raining. 

An example of a satellite derived precipitation area is shown in 
Figure 15. This is the corresponding area whose frequency distribution 
of brightness temperatures is shora in Figure 11; 14 June 1976 at 630 
GMT. Also indicated in Figure 15 are the rainfall amounts in hundredths 
of an inch recorded at the surface gages for the hour 600 to 700 GMT. 
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Figure 14^ Flow diagram of the decision process used to classify 
fields of view as raining or non-raining^ Tg is 
observed surface temperature, RHg is surface relative 
humidity, T is the mean of the Gaussian distribution 
of non-raining areas, and o is its standard deviation. 
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,e derived rain areas and amounts of rainfa 
ths of an inch) at surface observing rain- 
ir 14 June 1976 at 0630 GMT. Area of satel 
lervation is also indicated j 





4.0 RESULTS AI\’D DISCUSSION 


4.] Surface Rainfjaf'e Results 

Surface rain^’,aj'e obKerv'ations (sei- I’ir.uri: .10) have Ueen analv^-.i d 
for comparison with satellite derived prucii'ltation iri'(jmau-y. I’ij'urt;; 
If), 18, 20, 22, 24, 2b, 28 and 30 shov/ the results ol the surface pre- 
cipitation observations for June, July, Auj'ust and total for t'.e tiirce 
month period (day and night observations being kept separate) , Only 
surface observations which corresponded to the hour of a .satellite ol)- 

servation were tabulated. 

/ 

4.2 Results of the Satellite Derived Precipitation Frequency 

The results of the analysis described in Section 3.6 are shuv7n 
in Figures 17, 19, 21, 23, 25 and 27 and give the number of precipitation 
observations within 100 x 100 km squares for June, July and August, day 
and night observations being kept separate. Figures 29 and 31-give the 
precipitation observations as the percent of the total observations 
for the total three month period, again day and night observations 
are kept separate. 

4.3 Discussion 

4.3.1 The Climatology 

The relationship between the surface observed precipitation via 
hourly observations with raingages and satellite derived precipitation 
via microwave brightness temperature analysis is not straight forward. 

The tV7o systems are in fact measuring somewhat different quantities, 

Raingages measure the precipitation which reaches the surface at a 
single point leaving large areas unobserved. Also, hourly observations 
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Figure 20>. Same as Figure 16 but for July daytime orbits 























JULY NIGHT 



Figure 23^ Same as Figure 17 but for July nighttime orbits^ 
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AUGUST DAY 



Figure 24. Same as Figure 16 but for August daytime orbits 




Figure 25. Same as Figure 17 but for August daytime orbits„ 
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Figure 31o Sar-v* as Figure 29 but for the total three month 
nighttime orbitSo 
































give the integrated rainfall over a complete hour so that more than a 
single raining sy.stem may contribute to a gage’s recorded total. 

The satellite observed scattering by rain drops which may or may 
not roach the ground. It "sees" all areas, however it integrates c>very-- 
thing within one field of view (20 km x 42 km) into a single observation. 
Raining clouds must therefore be extensive enough or raining hard enough 
to affect the radiative condition of a large area and be observed by the 
satellite. The satellite also takes an effective "snapshot" of the 
area giving only an instantaneous view of the area’s precipitation. 

With the above differences in the measuring systems in mind, the 
major comparative features of the two analyses are discussed below. 

JUNE (10th through the 30th) 

Day - The absence of rain in Nebraska and Kansas are apparent in 
both the satellite and raingage observation. Relative maximums in the 
western Dakotas, central and southern Minnesota, northern Iowa and most 
of Illinois are also seen by both systems. 

Night - Minima seen in both the satellite and raingage observations 
are in the eastern Dakotas, Nebraska and Kansas, Illinois, Wisconsin 
and southern Missouri. Maxima are seen in the western Dakotas, north- 
eastern Minnesota, southern Iowa, and eastern Kansas. 

Day - Rain events observed by the gages are very scattered over the 
area. The satellite observed substantially more rain events than the 
gages and the spatial correlation between the two systems of observation 
is only generally good. The minima in the eastern Dakotas and central 
Minnesota are apparent in both systems. Maxima in central South Dakota, 


eastern Nebraska and northeastern Minnesota correlate well. The areas of 
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rain in northern North Dakota, eastern Kansas anJ Illinois as seen by 
the satellite are not apparent from surface raiagage observation. 

Night - The correlation of the midnight observations is much better 
in July than for the noon observations. Minima in central Minnesota, 

South Dakota, Missouri and Illinois correlate well between the two 
systems. Maxima in North Dakota, eastern Minnesota to northern Wisconsin, 
Iowa and western Kansas show up well in both systems as well. 

AUGUST 

Day - As in July the gage observations are highly scattered and the 
satellite observes more events than do the raingageso The spatial cor- 
relation between the systems, however, is better than for the July 
noon observations. Correlative minima are in northern Minnesota, 
northern Wisconsin, eastern South Dakota, Kansas and Missouri. Maxima 
in western North Dakota, eastern Nebraska, western Iowa and western 
Illinois are also highly correlated. 

Night - Good correlation is again apparent both spatially and in 
number of events for this nighttime case. Minima covering most of South 
Dakota, Kansas and Missouri are seen in both systems. The rain areas in 
North Dakota, Nebraska and northern Illinois are all seen in both the 
satellite observations and the raingages. The large number of rain areas 
as seen by the satellite in northern Iowa are not observed by the surface 
gages, however n 

THREE MONTH TOTAL 

Day - Spatial correlation is quite high although the satellite 
observes more rain events than do the gages. Figures 28 and 30 give 
rainfall in percent of total observations but due to July and August 
especially the number of events is quite high. The prominent minima 
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which show up are in eastern South Dakota, western Nebraska, most of 
Kansas, northern and eastern Missouri, eastern Iowa and northern Wis- 
consin. Maxima apparent in both systems are in western North Dakota, 
central South Dakota, central Nebraska, southeastern Missouri, western 
Iowa, northeastern and southern Minnesota, southern Wisconsin and most 
of Illinois. An extensive area of large percentages of rainfall is in 
northern North Dakota into Minnesota. This area contains the least 
dense raingage distribution (see Figure 10) so that a lot of scattered 
precipitation events could occur without being detected by any raingagesj 
Night - Spatial as well as qviantitative correlation between the 
two systems is very good. Minima are seen in eactern South Dakota, 
vjestern Minnesota, Kansas and Missouri all of which appear in both 
systems. The maxima which correlate well in central North and 
South Dakota, western Nebraska, southeastern Kansas, Iowa, eastern 
Minnesota and western Wisconsin. 

4.3.2 Regression of Satellite Observations on Raingage Observations 
Tabulation of the raingage observations into number of events 
in 100 X 100 km square areas as with the satellite data has been done 
to compare quantitatively the results from the two systems of observation. 
To examine how the satellite observed precipitation is related to the 
surface raingage observations a linear regression analysis has been 
done. 

The area has been divided into 4 quadrants and the results of the 
regression analysis in each quadrant are shown in Figures 32 and 33. 

The straight lines are the sample regression curves for each quadrant 
x^ith X representing the raingage observations and Y the satellite 
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y=NO. OF SATELLITE PRECIPITATION OBSERVATIONS 


Figure 32 ^ Results of regression analysis of the 4 quadrants 

for daytime cases. Straight lines are the regression 
curves and the hyperbolas indicate the 957 <, confiaence 
intervals of the regression curve. 
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The hyperbolas above and below the rej^reasion curves i 
the 95% confidence intervals of the sample repression curves. That is 
We are certain to 95% that the rep,ression curve lies vjithin the.su liiui.^ 
Also indicated on Figures 32 and 3 3 are the averapus of ttie two varialut 
(alvso indicated with an X on the repression curves), the sample* re- 
pression coefficient b, t^ which is a test parameter on the significance 
of b and p^ which is the probability of b actually being zero, Y is 
then independent of X. 

It is apparent that the best relationship between satellite and 
raingage observations is exhibited in the southern half of the area at 
both day and night (p^ is less than about 2% over the southern part 
of the area) . The northeastern quadrant shows by far the worst rela- 
tionship x\rith large p^ (40.6% night and >50% day). 

The southern half of the area is expected to be the best suited to 
the technique since x^armer surface temperatures provide a better back- 
ground v;ith xvhich to observe raining areas from the satellite. The 
colder surface temperatures to the north, particularly at night mask 
the raining areas. Also, the northeastern quadrant contains Lake 
Superior and Lake Michigan as Xv^ell ae numerous small lakes throughout 
Minnesota and Wisconsin. Although polarizing xc^ater surfaces have been 
eliminated from the analysis, small lakes which do not fill a field of 
viex \7 as x\^ell as rough wind blox\ai lakes can create ambiguities in clas- 
sification of raining areas. 

It can be said x^;ith confidence (within 2 to 3%) that indeed 
satellite observed rainfall is related positively to surface observed 
rainfall over the northx^est, southwest and southeast quadrants of the 
area during daytime cases and over the southern half of the area during 


nighttime cases. 
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4.3.3 Diurnal Variation 

VJallaco (1975) inventi^iated the diurnal nature oJ' preci])ltation 
frequency over the U.S. by examining, hourly precipitation objK‘rvatic)ii.s 
and has shown that over the north central U.S., including, much of the 
area of this study j. there is a pronounced near midnlpdit maximum in 
thunderstorm and heavy rainfall frequency during the summer season. 

The same diurnal character is observed in the hourly raingage observa- 
tions for June-August at least for midnight relative to noon. However, 
precipitation frequency derived from satellite microwave brightness tem- 
peratures shows more events around noon than near midnight. 

From the average values of the two variables in Figures 32 and 33, 
it is seen that for the nighttime case the raingages observe more events 
than the satellite and for the daytime case the satellite observes more 
events than the raingages. Therefore, it is believed that the satellite 
method is biased, observing more rain events than actually occur for 
the near noon case and observing fewer events for the near midnight case. 

The probable cause of the daytime bias is the larger variation "'n 
surface temperatures during the noon hours over the area which may be 
enhanced by cloud cover. That is areas overcast with non-raining clouds 
will have considerably lower surface temperatures than the surrounding 
sunlit areas causing their brightness temperatures to fall in the cold 
end of the frequency distribution and possibly below the threshold. 
Another possible cause may be that the satellite method during the day 
observes a large amount of rain which never reaches the ground but is 
evaporated into the warmer (relative to the midnight case) lower rela- 
tive humidity air below cloud base. Both of these causes include the 
existence of clouds which could also account for the maintenance of the 


spatial correlation between the two observin?; Hvtjte.’H'-. Jui'inr, the near 
noon casea . The "miaclasaif ied" (rain cioea not reach tiu- ground; event; 
by the satellite would tend to be near raiuinp, systems (other cloud:;) 
and simply expaird the areal extent of the "real" rain areas. 

The nighttime bias is due to the occurence of colder surface tem- 
peratures which mask the raining areas. From Figures 4 and 13, an area 
with a rainfall rate of 4 mm/hr has a brightness temperature of 255 K 
while a non-raining area with a surface temperature of 283 K (10"c) 
would also be expected to have a brightness temperature of 255 K. 

4.3j 4 Total Precipitatior 

Maps of total precipitation for June, July and August 1976 (Weekly 
Weather and Crop Bulletin) are shoTO in Figure 34. It is apparent that, 
as in the noon and midnight observations previously discussed, June was 
the rainiest month over the north central U.S. followed by July and then 
August. Furthermore, the major spatial features of the total precipita- 
tion are seen as well in the satellite derived precipitation. Maxima 
in central Iowa and northern North Dakota and Minnesota (especially in 
June and August) are apparent in both sets of figures. Also, minima 
throughout much of South Dakota, eastern Nebraska and Kansas as well 
as parts of Missouri and Wisconsin are generally correlated over the 
three month period. 







5.0 CONCLUSIONS 


lIicrov;ave ttnnperatures irom have* butai analy?:ud 

to clott'rmino raiulnj; a«jas over tho north central U.S. during, the montiKs 
viunu throuy.h Auj^uat 197b. The precipitation frequency of the area de- 
rived by the satellite has been compared witli that derived irom hourly 
surface rain^age observations. Both noon and midnight cases are seen 
to exhibit a positive relationship over much of the area between satellite 
observed rainfall events and surface raingage events. Only the northeast 
quadrant of the area shows little relationship. The satellite method 
is biased such that more events are observed during the daytime cases 
than actually occur at the surface and fewer events arc seen during the 
nighttime cases than are observed at the surface. It x^as also found by 
qualitative comparison that extreme maxima and minima in the total pre- 
cipitation of the area correspond to areas of maximum and minimum pre- 
cipitation frequency as derived by the satellite microwave analysis. 

It should be noted that although the thresholding technique des- 
cribed in Chapter 3 has x7orked Vv^ell in this study, general effectiveness 
of the technique should not be tacitly assumed for other areas at other 
times of the year. The technique works during the summer months because 
of hot surface brigheness temperatures which provide a background dis- 
tinguishable from the cold brightness temperatures of raining art^as^ 
Ambiguities occur when conditions change the radiative characteristics 
of the surface or intervening non-raining atmosphere such as to lox^rer 
the observed brightness temperature. For this reason supplemental infor- 
mation is needed to resolve as many ambiguities as possible; e.g. surface 
temperatures, cloudiness, dew formation, polarization, etc. The use of 
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man- interactive computer analyni.s ayiiCeniH will allow multi-parameter 
deciHi(m makinp, to be done efficiently. 

Precipitation frequency has beeit derived in this study, however 
other importaiit precipitation parameters may also be obtainable froi.i 
tills mode of analysis. Areal extant of rainfall events is observable 
via the satellite and is difficult at best to derive from synoptically 
spaced raingage measurements. Volumetric estimates of precipitation 
may be possible if rates can be assigned to the satellite derived 
rainfall areas, a difficult if not impossible task due to the low re- 
sponse in brightness temperature from a change in rainfall rate compared 
to changes in surface temperature and emittance , This may be done in- 
directly from comparison with surface raingage and radar observations. 
Studies of the distribution of rainfall for varying synoptic conditions, 
i.e. warm frontal, cold frontal, convective, etc., are also feasible 

Further theoretical and observational investigations are needed in 
the areas of microwave interaction with the atmosphere and how various 
surface types affect microwave emission. Studies utilizing extensive 
ground truth, bc.;h of precipitation and surface temperature, are neces- 
sary. Also, additlona] observing frequencies would largely decrease 
possible ambiguities. The Nimbus-7 satellite carries a Scanning Multi- 
channel Microwave Radiometer (SMMR)* which measures radiation at 6.6, 
10.7, 18, 21 and 37 GHz. Information expected from SMMR include sea 
ice parameters, soil wetness, snow accumulation over ice sheets as 
well as atmospheric total water vapor and non-precipitating liquid water 
to enhance rainfall estimations. An important advance in precipitation 




See Gloerson and Hardis (1978) o 


observation from satellites would be a j;eosynehronous orbiting, mierowave 
radiometer (Staeiin and Rosenkraiii',, 1978). This would allow lor 1 reciumii 
observation and thus, studies ol tiu\soscaU‘ pri ? ipi tatin^’; systems. 
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